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ABSTRACT 

We report an investigation of the extremely metal-poor and C-rich planetary nebula (PN) K648 in the globular 
cluster MIS using the UV to far-IR data obtained using the Subaru, HST, FUSE, Spitzer, and Herschel. We 
determined the nebular abundances of ten elements. The enhancement of E ([E/H]=h- 0.96) is comparable to that 
of the halo PN BoBnl. The central stellar abundances of seven elements are determined. The stellar C/O ratio 
is similar to the nebular C/O ratios from recombination line and from collisionally excited line (CEL) within 
error, and the stellar Ne/O ratio is also close to the nebular CEL Ne/O ratio. We found evidence of carbonaceous 
dust grains and molecules including Class B 6-9 pm and 11.3 ^m polycyclic aromatic hydrocarbons and the 
broad 11 /rm feature. The profiles of these bands are similar to those of the C-rich halo PNe H4-1 and BoBnl. 
Based on the theoretical model, we determined the physical conditions of the gas and dust and their masses, 
i.e., 0.048 Mq and 4.95 x 10“^ Mq, respectively. The observed chemical abundances and gas mass are in good 
agreement with an asymptotic giant branch nucleosynthesis model prediction for stars with an initial 1.25 Mq 
plus a 2.0 X10“^ Mq partial mixing zone (PMZ) and stars with an initial mass of 1.5 Mq without a PMZ. The 
core-mass of the central star is approximately 0.61-0.63 Mq. K648 is therefore likely to have evolved from a 
progenitor that experienced coalescence or tidal disruption during the early stages of evolution, and became a 
^1.25-1.5 Mq blue straggler. 

Keywords: ISM: planetary nebulae: individual (K648), ISM: abundances, ISM: dust, stars: Population II 


1. INTRODUCTION 

Planetary nebulae (PNe) represent a stage in the evolution 
of initial ^1-8 Mq stars. At the end of their evolution, such 
stars evolve into asymptotic giant branch (AGB) stars, then 
PNe, and finally white dwarves (WD). During this process of 
evolution, these stars eject a large fraction of their mass into 
the interstellar medium. The history of the progenitors is im¬ 
printed in the central star of the PN (CSPN) and the ejected 
gas. An investigation of the CSPN and the ejected material 
provides useful information to increase our understanding of 
stellar evolution, as well as the chemical evolution of galax¬ 
ies, i.e., how much of the mass of the star becomes a PN, 
which and how much of the elements are synthesized in the 
inner core of the progenitor, and how galaxies become chem¬ 
ically rich. The ejected gas in the PNe consists of both pro¬ 
cessed and unprocessed matter: primordial sources of proto¬ 
star cluster clouds or intracluster medium, pollution sources 
from highly evolved stars AGB and supernovae (SNe), and 
the result of stellar evolution processes (nucleosynthesized el¬ 
ements, molecules, and dust). Our understanding of the evo¬ 
lution of low-mass stars formed in the early Galaxy, as well 
as the chemical evolution of the Galaxy, can be enhanced by 
studying metal-poor PNe located in the Galactic halo. 

Eourteen Galactic halo PNe have be en identified 
since t he discovery o f K64 8 in M15 (e.g., [Howard et^ 
19971 iJacobvetall 119971 iPeauignot & Tsami^ 120051 
Pereira & Mirandall2007h . l^centlv. the number of detections 
has stea dily increased due to the Sloan Digital Sky Survey 
(SDSS) (I Yuan & Liull2013h . Eive PNe are located in the glob¬ 
ular clusters (GCs) M15 (K648), M22 (GJJCl and M2-29), 
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Pal6 (JaEul), and NGC6441 (JaFu2), and others are located 
in the Galactic halo field. The classification of PNe based on 
chemi cal abundances was originally proposed by iPeimberj 
(119781). and h as rec ently been revised and updated, e.g., 
lOuireza et al. (|20 07l). H alo PNe are classified as Type IV; 
specifically. iCosta et al.l (119961) indicated that halo PNe 
exhibit a large vertical distance from the Galactic plane ((z) = 
7.2 kpc) and large peculiar velocity relative to the rotation of 
the Galaxy ((AV) = 173 km s“', see their Table 6). Among 
halo PNe, H4- 1 (iTaiitsu & Otsukal 1201^ lOtsuka & Taiitsu 
2013), BoBnl (lOtsuka et al.ll201(11) 7 and K648 (iKwitter et ^ 
2003h are extremely metal-poor and C-rich (([Ar/H]) = 
-2.03, (C/O) = 14.49; this work); furthermore, there is an 
unresolved issue in terms of the chemical abundances; how 
did these progenitors evolve into C-rich PNe? The scientific 
backgr ounds of these PNe were explain ed by lOtsuka et al.l 
(1201(11) and by lOtsuka & Taiitsul (120131) . The progenitors 
of these three halo PNe were probably ^0.8 Mq stars, 
corresponding to the typical mass of turn-off stars in M15, 
because the [Ar/H] abundances as a metallicity indicator are 
si milar to the typic a l [Fe/H ] abundance in M15; according 
to iKobavashi et al.l (1201 ll) . [Ar/H]^-2.03 corresponds to 
[Fe/H]^-2.3. At least some of the stars of the Milky Way’s 
stellar halo were accreted along with their parent dwarf 
galaxies. BoBnl, a member of the oldest population in the 
Sagittarius dwarf spheroidal galaxy (IZiilstra et al.ll200^ and 
H4-1 in the halo field, might be younger than the classical 
Milky Way stellar halo population. 

For low-mass stars to evolve into C-rich PNe, a third 
dredge-up (TDU) is essential during the thermal pulse (TP) 
AGB phase. TDU conveys the He-shell reaction products, in¬ 
cluding C, O, Ne, and neutron («) capture elements, to the 
stellar surface. It is widely believed that >1-1.5 Mq stars ex- 
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Figure 1. The Ha image of K648 taken by the HST/WFFC2 with the F656N 
filter. North is up and east is to the left. The intensity is in erg s ~^ cm“^ 
Bright stars close to K648 are subtracted out. See Section [7^ for details of 
the methods employed. The locations of two RGB stars are indicated by the 
arrows. The reddening corrected magnitudes are 17.35 (B) and 16.78 (V) in 
RGBl and 17.54 (S) and 17.11 (V) in RGB2, respectively. 

perience TPU re.g.jLattanzT^ll987l : lKarakasl[2010l) . Recently 
iLugaro et alJ (1201 2^ reported the occurrence of TDU in ini¬ 
tial 0.9 Mq stars with a metallicity of Z = 10““*, although the 
minimum mass required for TDU depends on the model used. 
Even if TDU took place in the ~0.9 Mq progenitors, the post- 
AGB evolution of such low-mass stars toward the hot WDs is 
very slow. In addition, the ejected mass itself is very small, 
so it is difficult to observe them as visible PNe. Hence, the 
most likely explanation is that these progenitors gained mass 
via binary interactions to create new conditions for evolving 
into C-rich PNe. 

In view of the internal kinematics and nebular morphology, 
the progenitor of K648 appears to be a high-mass star. K648 
has bipolar and equatoria l outflows (iTaiitsu & Otsukall2006h 
and asymmetric nebulae (lAlves et al.l 120001k In Fig. [1] we 
show an Ha image of K648 obtained using the HSTfWFPC2. 
This image was processed using Lucy-Richardson deconvo¬ 
lution. K648 is composed of three parts: a very bright inner 
elliptical shell, an outer elliptical shell, and a bright arc on 
the northwestern limit of the major axis of the nebula, located 
just inside the edge of the outer bright elliptical shell. The 
arc is especially prominent in this object. A corresponding 
feature at the other end of the major axis does not appear to 
be present, although two fairly bright red giant branch (RGB) 
stars are unfortunately superposed at this location, making it 
difficult to resolve this feature. The locations of these RGB 
stars are indicated by the white arrows in the figure. The faint 
halo surrounding the outer elliptical shell extends to a radiu s 
of ~2.1" (not shown here, see Fig. 2 of I Alves et al.ll20()0h . 
The major axis of t he inner and outer shells is along the po¬ 
sition angle of -27°. iGarcfa-Segura et ^ (Il999h theoretically 
predicted that bipolar nebulae can be created in initial >1.3 
Mq single stars. We will explore the possibility of a binary 
system related mass-transfer activity suggested bv lAlves et alJ 
(I2000h . to solve the C abundance problem and the apparent 
contradiction in the evolutionary timescale. 

It would be interesting to study whether an increased mass 
star would evolve into a C-rich PN through such an evolu¬ 
tionary route. In AGB nucleosynthesis models, the predicted 


abundances, in particular «-capture elements, depend on the 
TDU efficiency, the number of thermal pulses, and the 
pocket mass. Any n-capture elements have not yet been de¬ 
tected in K648. The Ne ab undance is also sensitive to the 
amount of pocket mass (IShingles & Karak'a^l2013h . The 
Ne abundances can be easily determined using atomic gas 
phase emissions from the PNe rather than stellar absorption. 
To obtain a detailed view of the origin and evolution of K648 
through comparison with AGB nucleosynthesis models, we 
must accurately determine the abundances of C, O, Ne, and 
n-capture elements, and estimate the ejected mass. K648 is 
an ideal laboratory in which to investigate the evolution of 
low-mass metal-poor stars, as well as their nucleosynthesis. 
The reasons for this are first that the upper mass limit of stars 
in M15 is known (<1.6 Mq), and second that, because the 
distances are known with relatively little uncertainty, it is pos¬ 
sible to determine the core-mass of the PN as well as of the 
ejected mass. Study of K648 benefits not only understanding 
of the evolution of low-mass metal-poor stars, but also dust 
production in these stars. 

In this paper, we describe detailed spectroscopic analyses of 
K648 to investigate the origin and evolution of this PN based 
on an extensive set of spectroscopic/photometric data from the 
far-UV to far-infrared (FIR) regions of the electromagnetic 
spectrum. The remainder of the paper is organised as follows. 
In Section 2, we describe these observations using the Sub- 
aru/HDS, //5TAVFPC2/FOS/COS, ^pifzer/IRS/IRAC/MIPS, 
and Herschel/PACS. In Section 3, we provide the elemental 
abundances of the nebula and the CSPN, as well as the physi¬ 
cal properties of the CSPN. We determined the abundances of 
the 10 elements of the nebula of K648, including the first mea¬ 
surements of the n-capture element fluorine ( F) in this PN. Us¬ 
ing th e spectrum synthesis code TFUSTY (iFanz & Hubenvl 
HqqI, we determined the abundances of 7 elements of the 
CSPN and the core-mass of the CSPN. We also report the 
C-rich dust features found in the SpitzerflRS spectrum. We 
constructed a self-consistent model, whereby the predicted 
spectral energy distribution (SED) fits the observations and 
accordingly estimated the mass of eject ed gas and dust usin g 
the radiative transfer code CFOUDY (iFeriand et al.lll998ll . 
In Section 4, we compare the elemental abundances of K648 
with those of H4-1 and BoBn I. We discuss the origin and evo¬ 
lution of K648 by comparing the predicted elemental abun- 
dances, the final cor e-mass, and the ejected mass reported by 
iFugaro et al.l (120121) with our determined values. A summary 
is presented in Section 5. 

2. OBSERVATIONS & DATA REDUCTION 
2.1. HDS observations 

Optical high-dispersion spectra we re obtained using a 
High-Dispersion Spectrograph (HDS; iNoguchi et al.l 1200^ 
attached to the Nasmyth focus of the 8.2-m Subaru telescope 
on 2012 June 28 (Program ID: S12A-078, PI: M. Otsuka). 

The weather conditions were stable and clear throughout 
the night, and the seeing was ^0.5" measured using the guider 
CCD. An atmospheric dispersion corrector (ADC) was used 
to minimize the differential atmospheric dispersion through¬ 
out the broad wavelength region. The slit width was set to 
1.2" and the slit length was set to 6"; these settings allowed us 
to reduce contamination from nearby stars. We selected 2x2 
on-chip binning. The resolving power {R) was 33 500, deter¬ 
mined from the average full-width at half-maximum (FWHM) 
of over 600 Th-Ar comparison lines. Blue-cross and the red- 
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Wavelength (rest) (A) 

Figure 2. The de-reddened HDS spectrum of K648. The wavelength is cor¬ 
rected to the rest wavelength in air. 


cross dispersers were employed to obtain the 3620-5400 A 
spectrum (blue spectrum) and the 4320-7140 A spectrum (red 
spectrum), respectively. We set the position angle to -27° us¬ 
ing an image de-rotator. The total exposure times were 7200 
s for the blue spectrum and 9000 s for the red spectrum, re¬ 
spectively. Flux calibration, blaze function correction and air- 
mass correction were carried out by observing the standard 
star BD-i-28° 4211 twice at different airmasses for each blue 
and red spectrum. 

Data reduction was carried out using the Echelle Spectra 
Reduction Package ECHELLE in IRAF0. We generated a 
single 3620-7140 A spectrum by combining the blue and the 
red spectra after scaling the flux density of the blue spectrum 
by a factor of 1.06 to match that of the red spectrum. The re¬ 
sulting signal-to-noise (S/N) ratio was >50-130 for the con¬ 
tinuum of this single 3620-7140 A spectrum. 

Eigure |2] shows the resulting spectrum, which is corrected 
for interstellar extinction (see the following section). The ob¬ 
served wavelength was corrected to the averaged line-of-sight 
heliocentric radial velocity of-116.89±0.41 km s“' (the root 
mean square (RMS) of the residuals was 4.15 km s“') among 
all lines detected in the HDS spectrum (122 lines). 


2 . 2 . Interstellar reddening correction of the HDS spectrum 

The line-fluxes were de-reddened using the follow expres¬ 
sion: 

/(A) = F{\) X io°<h/3)(i+/(A))^ ( 1 ) 

where /(A) is the de-reddened line flux, F{X) is the observed 
line flux, /(A) is the interstellar extincti on function at A com- 
puted by the reddening law reported bv ICardelli et al.l (119891) 
with Rv = 3.1, and c(H/3) is the reddening coefficient at H/3. 
Our measured E(H^) was 1.70x 10“^^±3.14x erg s“' 
cm“^ in the HDS spectrum. Hereafter, X(-Y) corresponds to 

^ IRAF is distributed by the National Optical Astronomy Observato- 
lies, operated by the Association of Universities for Research in Astron¬ 
omy (AURA), Inc., under a cooperative agreement with the National Science 
Foundation. 


Xxl0“'^. We computed c(H/3) by comparing the observed 
B aimer line ratio s of H 7 and Ha to H/3 wi th the theoretical 
ratios reported by I Storey & Hummeil (119951) with an electron 
temperature of T^= 10"^ K and an electron density of n^ =10"^ 
cm“^ with the assumptions of Case B. The values of c(H/3) 
were 0.121±0.027 from the E(H 7 )/E(H^) and 0.148±0.008 
from the F{Ha)/F{Hf]) ratios. We used an average value of 
c(H/3) = 0.135 ±0.017 for the interstellar reddening correc¬ 
tion. 

2.3. Emission-line flux measurements with the HDS spectrum 

The detected emission-lines are given in the Appendix (see 
Table ©. Eor the flux measurements, we applied multiple 
Gaussian component fitting. We list the observed wavelength 
and de-reddened relative fluxes for each Gaussian component 
(indicated by Comp.ID number in the fourth and the eleventh 
columns of Table lAl in the Appendix), with respect to the de- 
reddened flux of 100. /(A) for each wavelength is also 
listed. Most of the line-profiles of the detected lines can be 
fitted using a single Gaussian component. Eor the lines com¬ 
posed of multiple components, e.g., [Oil] A3726.03 A, we 
list the de-reddened relative fluxes of each component, as well 
as the sum of these components (indicated by Tot.). 

We supplemen t ed ou r HDS data with the data given by 
iTaiitsu & Otsukal (l2006h to calculate the Ar^^ abundance us¬ 
ing [Arm] A7135 A, ^(.([Oll]) and «e([Oll]) by combining 
[On] A A 7320/30 A with [Oil] AA 3726/29 A, and 7;(Hei) 
using He IA 7281 A. 

2.4. HST/WFPC2 photometry and the Ha/H(3 fluxes 

In the FOS UV-spectrum (see the following section), no 
H I or He II nebular lines are required to normalize the C III] 
and [C II] fluxes to the H/3 flux. Therefore, we measured 
the H/3 flux of the entire nebula and scaled the FOS flux 
density to tune the UV flux densities at bands including the 
C III] AA 1906/09 A and the [C ii] A 2323 A lines. The H/3 flux 
of the entire nebula is also necessary to normalize the fluxes 
of the lines detected in the Spitzer/IRS spectrum (see the fol¬ 
lowing section). Broadband fluxes are required to estimate 
the core-mass of the CSPN and to constrain the incident SED 
of the CSPN and the emergent spectra predicted by the neb¬ 
ular model. Eor this purpose, we used HSTPHide Eield and 
Planetary Camera 2 (WPPC2) photometry using eight broad¬ 
band and P656N Alters, which are available in the Mikulski 
Archive for Space Telescopes (MAST). 

We reduced the WPPC2 data (IDs: 10524 and 11975, 
PEP. R. Prancesco; ID:6751, PI: H. E. Bond) using the stan¬ 
dard HST pipeline and MultiDrizzle on PYRAP to re¬ 
move cosmic-rays and improve the angular resolution. Pirst, 
we removed nearby stars using empirical point-spread func¬ 
tions (PSPs) generated from IRAP/DAOPHOT. We then 
measured the count rates (cts) within an aperture radius of 
2.1". We defined the background sky as being represented by 
an annulus centered on the CSPN with inner radius of 3.2" 
and outer radius of 4.2". Pinally, we converted the cts into 
the flux densities using the PHOTPLAM values in erg s“* 
cm“^ cts“'. The resulting flux densities F\ and the cor¬ 
responding de-reddened data I\ are listed in the fourth and 
fifth columns of Tabled] 

To measure the Ha flux of the entire nebula using the 
P656N flux density, it is necessary to remove the contribu¬ 
tions of both the local continuum and the [N ll] A 6548 A 
line. However, this procedure entails a number of problems 
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Table 1 

//5TAVFPC2 photometry data for K648. 


Filter 

-^cen. 

(A) 

AA 

(A) 

CSPN+Nebula 

CSPN 

Prop.ID 

Fa 

(erg cm“^ A“^) 

A ^ 

(erg cm“^ A“^) 

Fa 

(erg cm“^ A“^) 

^A 

(erg cni“^ A“^) 

F160BW 

1515.16 

188.43 

1.05(-13)±2.54(-15) 

2.09(-13)±5.09(-15) 



11975 

F170W 

1820.78 

285.52 

9.72(-14)±9.01(-16) 

1.89(-13)±1.75(-15) 

9.54(-14)±7.67(-15) 

1.85(-13)±1.49(-14) 

11975 

F255W 

2598.57 

171.21 

3.01(-14)±3.95(-16) 

5.29(-14)±6.95(-16) 



10524 

F300W 

2989.04 

324.60 

2.19(-14)±2.16(-15) 

3.54(-14)±3.48(-15) 

1.43(-14)±1.28(-15) 

2.30(-14)±2.06{-15) 

11975 

F336W 

3359.48 

204.49 

2.32(-14)±3.30(-15) 

3.56(-14)±5.06(-15) 

1.76(-14)±4.94(-16) 

2.71(-14)±7.58(-16) 

6751 

F439W 

4312.09 

202.32 

1.11(-14)±4.19(-15) 

1.58(-14)±5.98(-15) 

9.97(-15)±3.46(-16) 

1.42(-14)±4.93{-16) 

6751 

F547M 

5483.88 

205.52 

4.62(-15)±9.57(-16) 

6.01{-15)±1.24(-15) 

3.84(-15)±1.04(-16) 

4.99(-15)±1.35(-16) 

6751 

F814W 

7995.94 

646.13 

1.57(-15)±2.97(-16) 

1.84(-15)±3.47(-16) 

1.24(-15)±5.01(-17) 

1.45(-15)±5.86(-17) 

6751 

F656N 

6563.76 

53.78 

1.04(-13)±4.33(-16) 

1.28{-13)±5.37(-16) 



6751 


Note. — Fx and lx are the reddened and de-reddened flux densities, respectively. We used the reddening law reported bv ICardelli et alj 119891) for interstellar reddening correction 
with /?v = 3.1 and E(B -V) = 0.092. 


(see e.g.. lLuridiana et aP (120031) for a thorough discussion of 
the pitfalls and uncertainties in determining line-fluxes from 
HST/WFPC2 images for the PN NGC6543). An alternative 
is to use an intermediate step of computing an equivalent Ha 
flux, that is, using the Spitzer HI Pfa and Hua recombina¬ 
tion lines. This method also has problems of contamination 
due to the 7.7 pm PAH feature, as well as the broad 11 pm 
feature, and the [Nell] A 12.80pm line in the Spitzer spec¬ 
tra. Rather than employing the above HST Ha flux extraction 
method, or the intermediate step of using Spitzer HI lines, 
we used the HSTPNFPCl F656N band flux density itself, 
i.e., Fx{HST,¥656N), which includes the Ha flux, the local 
continuum, and the [N ll] A6548 A line within the F656N fil¬ 
ter band, together with our Subaru/HDS spectrum. The ad¬ 
vantage of this approach is that it is possible to extract the 
Ha flux without contamination from nebular and stellar con¬ 
tinuum and [N ll] A6548 A. This method was applied in our 
previous work on the PN Ml-11 (lOtsuka et al.ll2013l) . Tak¬ 
ing into account the F656N filter transmission characteristics, 
we compared the F\(HST,F656N) with the counterpart Sub¬ 
aru/HDS scan spectrum, i.e., F\(HDS,F656N). The scaling 
factor /A(//57’,F656N)/F\(HDS,F656N) = 1.428 was deter¬ 
mined, and was applied to the Subaru/HDS spectral line fluxes 
to analyze both the spectra on an equal footing. After applying 
the scaling factor, the HDS fluxes should be F’(Ha) = 2.42(- 
12)±4.44(-14) erg s”' cm'^ and F{n(3) = 7.83(-13)±1.05(- 
15) erg s“* cm“^. A simple comparison of these scaled HDS 
data with the measured HST data shows very small devia¬ 
tions, i.e., 0.18%, &. 0.13%, corresponding to the Ha and H/3 
fluxes measured from non-scaled HDS spectrum. The un¬ 
certainties of our measurements are much s maller than the 
estima ted uncertainty of ^^10% reported bv iLuridiana et all 
(l2003h . These reduced errors may be coincidental and the ac¬ 
tual errors could be larger than our estimation; however, the 
errors in ou r analysis appea r to be smaller than the estimates 
reported bv ILuridiana et al.l (l2003h . The scaling factors also 
give ratios of Pfa and Hua with the above H/3 fluxes that are 
consistent with the theoretical values (see Section 2.7). Note 
that the SpitzerllF.S> spectra were obtained using a wider slit 
width, which was sufficient to cover the entire K648 nebula. 

2.5. HST/FOS UV-spectrum 

To calculate the and C*' abundances using the 

C III] A A 1906/09 A and the [C ll] A 2323 A lines, we analyzed 
the archival HST/FOS spectrum (The Faint Object Spec¬ 
trograph), which was obtained on 1993 Nov 18 (Prop.ID; 
3196, PI: H. Ford), and was downloaded from MAST. We 


Table 2 

The detected lines in the HST/FOS spectra. 


•^obs 

Ion 

“^lab 

Comp. 

/(A) 

i{\r 

(A) 


(A) 



[/(H/3) = 100] 

1906.83 

Cm] 

1906/09 

1 

1.256 

334.984±17.038 

2326.45 

[Cii] 

2323 

1 

1.392 

17.091±1.634 


^ We usedFtH/?) = 7.83(-13) ih 1.04(-15) erg cm““, which was measured based on 
the HSTfWFPCl F656N image and the observed F(Ha)//'(H/3) ratio of 3.097 measured 
using HDS spectra. 


used the data sets Y1C40103P, Y1C40104T, Y1C40105T, and 
Y1C40106T. 

We scaled the flux density to fit the F160BW, F170W, and 
F255W bands listed in the fourth column of Table [T]using the 
relevant transmission curves (scaling factor = 0.837). Using 
the F’(H^) = 7.83(-13)±l.p5(-15) erg s'^ cm'^, we normal¬ 
ized the C III] AA 1906/09 A and the [C ll] A 2323 A fluxes. 

2.6. FUSE and HST/COS UV-spectra 

We analyzed archival UV spectra of K648 from MAST 
to calculate the elemental abundances in the photosphere of 
the CSPN and determine the parameters required to calcu¬ 
late the stellar radius, surface gravity logg, effective tem¬ 
perature Tgff, and the current core-mass of the CSPN. The 
920-1180 A and the 1170-1780 A spectra were obtained us¬ 
ing the Far Ultraviolet Spectroscopic Explorer (FUSE) on 
2004 Nov 1 (data set: D1570101000, PI: Dixon) for and the 
HST/Cosmic Origins Spectrograph (COS) on 2013 Nov 13 
(data set; LB2402010/20; Prop-ID:l 1527, PI; J. Green). We 
generated the FUSE, HST/COS, and HDS spectra normalized 
to the flux density at a continuum of 1.0 using IRAF/S PLOT. 

2.1. Spitzer/IRS mid-infrared spectra 

We reduced the a rchive data obtaine d using the Infrared 
Spectrograph (IRS; iHouck et al.l 120041) with the SL (5.2- 
14.5 pm and a slit dimension of 3.6"x57"), SH (9.9-19.6pm, 
4.7"xll.3"), and LH( 18.7-37.2pm, ll.l"x22.3") modules 
(AOR Keys: 15733760 for the SL and 18627840 for the SH 
and LH spectra; Pis; R. Gehrz and J. Bernard-Salas, respec- 
tively). We used the data reduction packages SMART v.8.2.5 
(iHigdon et al.l2004h and IRSCLEAN provided by the Spitzer 
Science Center. Lor the SH and the LH spectra, we subtracted 
the background sky using the offset spectra. We scaled the 
flux density of the SL data to that of the SH &. LH data in the 
overlapping wavelength region. The remaining spikes in the 
spectra were removed manually. 
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Wavelength (/i.m) 

Figure 3. The Spitzer/IRS spectrum of K648. The detected gas emission 
lines listed in Tableware also indicated here. The 6-9 /rm PAH band, the 
11.3 fim PAH emission, and the broad 11 pm band are also indicated. The 
line-profile of the [S III] A 18.71 pm is zoomed in the inner box. The position 
of this line in the laboratory is indicated by the vertical red line. 


Table 3 

The detected atomic lines in the Spitzer spectra. 


-''lab 

(^m) 

Ion 

/(A) 

F{\) 

(erg cm“^) 

/(A) 

[/(H/3) = 100] 

lAl 

Hi 

-0.990 

3.23(-14)±1.56(-16) 

3.15±0.15 

8.99 

[Ar III] 

-0.959 

3.29(-15)±4.84(-16) 

0.32±0.05 

10.51 

[SIV] 

-0.959 

1.08(-14)±4.46(-16) 

1.06±0.07 

12.37 

Hi 

-0.980 

1.04(-14)±5.17(-16) 

1.02±0.07 

12.80 

[Ne II] 

-0.983 

1.53(-13)±1.08(-14) 

14.98±1.28 

15.55 

[Nelll] 

-0.985 

1.18(-13)±1.42(-14) 

11.54±1.49 

18.71 

[Sill] 

-0.981 

1.36(-14)±1.98(-15) 

1.33±0.20 

33.47 

[Sill] 

-0.993 

6.27(-15)±1.12(-15) 

0.61±0.11 


The resulting spectrum is shown in Fig. [3 iBover et'dl 
(I2006^ reported the SL spectrum for K648 only. Therefore, 
the spectrum at longer wavelengths (i.e., beyond ^^14.5 ^m) 
is shown here for the first time. The line-profile of the 
[S III] A 18.71 /im, which is faint in K648 and also a impor¬ 
tant diagnostic line, is also shown in the inner box. 

The line fluxes of the detected atomic lines are listed in 
Table [3] We corrected for the interstellar reddening using 
Equation (1 ) and t he interstellar extinction function given by 
iFliiks et al.l (QHl. We computed c(H/3) = 0.12±0.02 by 
comparing t he theoretical /(H 1 7 47 rtm )//(H/3) ratio of 3.15(- 
2) given by lStorev & Humni^ (Il995h for T^ = 10"^ K and 
= 10"* cm“^ under the assumptions of Case B. He re, we used 
F{\{j3) = 7.83(-13) erg s“* cm“^ (see Section lZ^ . This result 
appears appropriate as the measured /(Hi 12.37pm)//(H/3) 
= 1.02(-2 ) is in good agreement with the theoretical data 
(1.051-2'). IStorev & Hummeill 19951) . 

K648 exhibits the 6-9 pm polycyclic amorphous carbon 
(PAH) band and the broad 11 pm feature. These two features 
are frequently seen in C-ric h PN e. We will discuss the details 
on these features in Section lT^ 

2.8. Spitzer/IRAC/MIPSphotometry 

To provide a constraint in the SED fitting at mid-infrared 
(MIR) wavelengths, we reduced archival Spitzer MIR im- 
ages obtained u sing the Infrared Array Camera (IRAC; 
lEazio et a l]|2004|) and the Multiband Imaging Spectrometer 
(MIPS: iRieke et ^l2004ft . We downloaded the basic cali¬ 
brated data and reduced it using MOPEX, which is provided 


Table 4 

Spitzer/lRAC/MIPS and Herschel/PACS photometry of K648. 


Band 

-^cen. 

(/^m) 

AA 

(Aim) 

Fx 

(erg cm“^ 

AORKEY(Spitzer)l 

OBSlD{Herschel) 

IRAC-chl 

3.51 

0.68 

1.24(-12)±1.53(-13) 

12030208 

IRAC-ch2 

4.50 

0.86 

6.16(-13)±4.80(-14) 

12030208 

IRAC-ch3 

5.63 

1.26 

4.78(-13)±4.51(-14) 

12030208 

IRAC-ch4 

7.59 

2.53 

4.43(-13)±2.82(-14) 

12030208 

MIPS-chl 

23.21 

5.30 

5.95(-14)±1.44(-15) 

12030464 

PACS-B 

68.92 

21.41 

1.86(-15)±3.82(-17) 

1342246710/11/12 

PACS-R 

153.94 

69.76 

3.40(-16)±4.00(-17) 

1342246710/11/12 


by the Spitzer Science Center, to obtain single mosaic images 
for each band. 

We carried out PSE fitting photometry of the IRAC images 
using IRAE/DAOPHOT. We adopted the position of K648 
measured in the //57'/F656N image and corrected the flux 
densities measured using PSE photometry by aperture pho¬ 
tometry of the PSE stars. Eor the MIPS 24 pm image, we 
measured the total count within a 1" radius region, and sub¬ 
tracted the background represented by the annulus centered 
on the PN with 20" inner and 38" outer radii, respectively. 
We used an aperture correction factor of 1.61, as listed in the 
MIPS instrument hand book. The measured fluxes are listed 
in Table |4] 

2.9. Herschel/PACS photometry 

By combining the MIR data from the Spitzer and the FIR 
data from the Herschel, we attempted to trace the ejected 
mass of K648 during the last TP as accurately as possi¬ 
ble. For this purpose, we analyzed archived 70 pm (PACS- 
B) and 160 pm images (PACS-R) obtained using the Pler- 
.ic/rel/Photodetecting Array Camera and Spectrometer (PACS; 
iPoglitsch et al.ll2()T^ . 

We downloaded the reduced PACS data of K648 (OBSID: 
1342246710/11/12, PI: M. Boyer) from the Herschel Science 
Archive (HSA). The PACS images are shown in Fig. |4] For 
comparison, we also show the //5'7'/F656N images. The plate 
scale of the image shown in Fig. HJa) is 0.025" pixel"' and 
that of the //5'7'/F656N image shown in Fig.Ub) corresponds 
to that of PACS-B. Fig.lUc) shows PACS-B data at 1" pixel"', 
and Fig.lUd) shows the plate scale of PACS-R at 2" pixel"*. 
The most likely position of K648 was determined using the 
//57’/F656N image, and is indicated by the white crosses. The 
light from K648 is partially contaminated by nearby stars. 

We used IRAE/DAOPHOT to measure the flux densities 
within a radius of 2 pixels in both the PACS-B and PACS- 
R bands. We regarded the median count within the annulus 
centered on the PN with an inner radius of four pixels and an 
outer radius of five pixels as the background. We corrected 
the measured flux densities of K648 using aperture photom¬ 
etry with correction factors of 4.29 for PACS-B and 3.83 for 
PACS-B0. The measured flux densities are summarized in Ta¬ 
ble® 

3. RESULTS 

3.1. Emission-line analysis 
3.1.1. CEL diagnostics 

® These coirection factors were computed using the table of encircled en- 
ergy fractions as a function of the radius of the aperture for the PACS filter 
bands in the NASA Herschel Science Center. 
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Figure 4. The HSTf¥65(i and the Herschelf^AC-S images. The size of each panel is 17" x 17". North is up and east is to the left. The plate scale of the image in 
(a) is 0.025" pixel“^. The plate scale of the HST/F656N image in (b) corresponds to that of PACS-B (1" pixel”*). K648 is located in the center of each image. 
In the PACS-B and PACS-R images, the location of K648 is indicated by the crosses. 



Figure 5. An rie-Te diagram. Each curve is labeled with an ID number given 
in Table[5] The solid lines indicate diagnostic lines of Te, whereas the broken 
lines indicate diagnostic lines of n^. 


Table 5 

Plasma diagnostics. 


ID 

Diagnostic 

Value 

He (cm 

(1) 

(2) 

(3) 

(4) 

(5) 

[S ll](A6716)/(A6731) 

[0 ll](A3726)/(A3729) 

[S III](A18.7 fim)/(X33.5 Atm) 

[0 u](A3726/29)/(A7320/30) 

[Cl lll](A5517)/(A5537) 

0.658±0.023 

1.852±0.076 

2.178±0.523 

8.576±0.108“ 

0.749±0.120 

25304330 

34304470 

511042100 

78904130 

713043170 


Balmer decrement 


7500-10000 

ID 

Diagnostic 

Value 

T,(K) 

(6) 

(7) 

(8) 
(9) 

[O III](A4959+A5007)/(A4363) 

[Ne m](A15.5/tm)/(A3869+A3967) 
[N ll](A6548+A6583)/(A5755) 
[Arm](A8.99Atm)/(A7135) 

108.649±4.425 

0.882±0.115 

79.200±9.333 

0.844±0.125 

12 3504190 
11 0904450 
10 3804530 
102704900 


He I(A5876)/(A4471) 

He l(A6678)/{A4471) 

He l(A7281)/(A5876) 

He l(A7281)/(A6678) 

3.01940.033 

0.83740.012 

0.04040.001 

0.14540.003 

42704300 

71004760 

63604150 

66804130 


(Balmer Jump)/(H11) 

U. lU2±U.UUt) 

11 630±93U 


^ The recombination contribution is corrected for the [OII] A A 7320/30 A lines. 


In the following analysis using CELs and RLs, we used 
the transition probabilities, collisional impacts, and recombi¬ 
nation coefficients listed in Tables 7 and 11 of iQtsuka et al.l 

dloToh . 

The electron temperatures and densities were determined 
using a variety of line diagnostic ratios by calculating the state 
populations using a multilevel atomic model. The observed 


diagnostic line ratios are listed in Table |5] where the num¬ 
bers in the first column indicate the ID of each curve in the 
ne-Te diagram shown in Fig.|5] The second, third, and final 
columns in Table |5] show the diagnostic lines, line ratios, and 
the resulting and Te, respectively. We obtained nine diag¬ 
nostic line ratios with different ionization potentials (IPs) in 
the range 10.4 eV ([S ll]) to 41 eV ([Ne III]), and determined 
a suitable and combination for each ion. 

For the [O ll] AA 7320/30 A lines, we eliminated the recom¬ 
bination contamination due to O^ '*' using the following expres¬ 
sion, which is given by iLiu et al.l (120001) : 

/fi([OII]AA7320/30) 

/(H/3) ^'^^(.lOV 

Using the ionic abundances derived from the recombina¬ 
tion O II A 4641.8 A line and with T). = 11 650 K, based on 
the Balmer jump discontinuity (see the following section), we 
found that /r([Oii] AA 7320/30) = 0.12±0.02. As we could 
not detect the N II and the pure O III recombination lines, 
we were unable to estimate the contribution of N^'*' to the 
[N II] A 5755 A line nor that of to the [O III] A 4363 A line. 

First, we computed with = 10000 K for all 
density diagnostic lines. ^..([Ne III]), ^^([Olll]), and 

^^([Arlll]) were calculated using = 6100 cm“^, which 
is the averaged between «e([Cllll]) and ^(.([Slll]). 
We calculated ^^([N ll]) using the «e([Oll]) determined 
from the [Oil] /(A3726)//(A3729) ratio. We used [Oil] 
/(AA3726/29)//(AA7320/30) as a density indicator for the 
^4500 cm“^ region, which is larger than the critical density 
of [O II] A 3726 A. 

Our values of 7) a nd iie are comparable to those reported by 
iKwitter et all (l2()0l . i.e., 7;([0 III]) = 11 800 K, 7;([N ll]) = 
9200 K, and n,([S ll]) = 1000 cnr^. 

3.1.2. RL diagnostics 

We calculated 7]; using the ratio of the Balmer dis continuity 
to /(H ll). We employed the method reported bv iLiu et all 
(120011) to calculate the electron temperature re(BJ). 

We calculated the He I electron temperatures using the four 
different Tp /He l) line ratios and th e emissivities of these He I 
lines from iBeniamin etTP (Il999h . in the case of = 10^ 
cm“^. 

The intensity ratio of a high-order Balmer line H« (where n 
is the principal quantum number of the upper level) to a lower- 
order Balmer line is also sensitive to the electron density. The 
ratios of higher-order Balmer lines to H/3 are plotted in Fig. 
along with theoretical values from IStorev & Hummed (119951) 
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Figure 6. The intensity ratio of the higher-order Balmer lines to (Case 
B assumption). The theoretical intensity ratios (dotted curve and triangles) 
are given for Te = 11 650 K (determined from the Balmer Jump) and rie = lO'^ 
cm“^. 



Figure?. The 4770-4890A Subaru/HDS spectrum of K648. The wave¬ 
length is connected to the rest wavelength in air. The locations of the 
[FII] AA 4790/4869 A lines in are indicated by the vertical red lines. The 
position of the [Felll] A4881 A is indicated by the vertical black line. The 
nan'ow absorptions are from cool stars, possibly two nearby RGB stars (see 

Fig.C). 


Table 6 

The adopted electron temperatures and densities. 


71 (K) 

(cm ^) 

Ions 

10380 

2530 

s+ 

10380 

3430 

C+,N+,0+(3726,29 A),E+,Fe2+ 

10380 

7840 

0+(7320,30A) 

10270 

6100 

C2+,Ne+,s2+,CF+,AF+ 

11090 

6100 

Ne2+ 

12 350 

6100 

02+,s3+ 


for and = 10000 cm“^. We ran small-grid calcu¬ 

lations to determine «£ in the range 5000-12500 cm“^, and 
found that the models in the range of = 7500-10000 cm“^ 
provided the best fit to the observed data. Te and deter¬ 
mined using the RL diagnostics are summarized in Table |5] 

3.1.3. CEL ionic abundances 


Table 7 

Ionic Abundances from CELs. 


x™+ 

Alab 

7{Alab) 

Xm+/H+ 

c+ 

2323 A 

1.64(+1)±1.63(0) 

2.55(-5)±7.28(-6) 

C2+ 

1906/09 A 

3.35(+2)±1.70(+l) 

6.91M)±3.23(-4) 

N+ 

5754.64 A 

4.31(-2)±2.44(-3) 

4.67(-7)±1.02(-7) 


6548.04 A 

9.01(-l)±1.23(-2) 

4.93(-7)±5.89(-8) 


6583.46 A 

3.18(0)±4.10(-2) 

5.89{-7)±7.02(-8) 

5.68(-7)±6.77(-8) 

0+ 

3726.03 A 

1.74(+1)±2.99(-1) 

1.34{-5)±2.50(-6) 


3728.81 A 

9.39(0)±3.50(-l) 

1.34(-^5)±2.60(-6) 


7320/30 A 

3.12(0)±3.91(-2)“ 

1.79C5)±4.40(-6) 

1.34(-5)±2.54(-6) 

0^+ 

4363.21 A 

2.78(0)±2.56(-2) 

4.03{-5)±3.29(-6) 


4931.23 A 

3.84{-2)±4.83(-3) 

5.12{-5)±6.80(-6) 


4958.91 A 

7.50(+l)±7.37(0) 

3.90{-5)±4.18(-6) 


5006.84 A 

2.27(+2)±9.46(0) 

4.09(-^5)±2.45(-6) 

4.05(-5)±2.89(-6) 

F+ 

4789.45 A 

1.10(-l)±3.67(-3) 

6.67(-8)±1.02(-8) 


4868.99 A 

2.96(-2)±3.36(-3) 

5.75{-8)±1.08(-8) 

6.47(-8)±1.03(-8) 

Ne+ 

12.80 fim 

1.50(+1)±1.28{0) 

2.01(-5)±1.95(-6) 

Ne2+ 

3869.06 A 

9.94(0)±1.24(-1) 

7.28{-6)±1.00(-6) 


3967.79 A 

3.15(0)±4.35(-2) 

7.65(^)±1.06(-6) 


15.55 /rm 

1.13(+1)±1.36(0) 

7.48(^)±9.79(-7) 

7.42(-6)±9.98(-7) 

S+ 

6716.44 A 

8.73(-2)±2.23(-3) 

6.72(-9)±7.79{-10) 


6730.81 A 

1.33(-l)±3.21(-3) 

6.73(-9)±7.48(-10) 

6.72(-9)±7.60(-10) 

S2+ 

6313.1 A 

1.19(-l)±5.12(-3) 

2.52(-7)±7.30(-8) 


18.71 fim 

1.33(0)±2.04(-1) 

2.10(-7)±3.52(-8) 


33.47 fim 

6.12(-1)±1.13(-1) 

2.10(-7)±4.16(-8) 

2.12(-7)±3.93(-8) 


10.51 fim 

1.06(0)±6.66(-2) 

3.35(-8)±2.11(-9) 

cF+ 

5517 . 72 A 

2.12(-2)±2.71(-3) 

3.15(-9)±7.86(-10) 


5537.89 A 

2.83(-2)±2.73(-3) 

3.17(-9)±7.38(-10) 

3.16(-9)±7.59(-10) 

Ai2+ 

7135.79 A 

3.84{-l)±6.20(-3) 

3.32(-8)±5.73(-9) 


8.99 Aim 

3.24(-l)±4.77(-2) 

3.41(^8)±5.32(-9) 

3.36(-8)±5.54(-9) 

Ee2+ 

4701.53 A 

2.22{-2)±3.11(-3) 

2.37(-8)±4.80(-9) 


4881.11 A 

5.09{-2)±3.55(-3) 

2.75{-8)±4.59(-9) 

2.63(-8)±4.65(-9) 


^ Corrected recombination contribution for the [OII] AA 7320/30 A lines. 


We obtained the following 14 ionic abundances: N'*', 

0+.2+^ P+, Ne+’2+, S+^2+.3+^ cj2+^ ^j.2+ pg2+ ppg 

dances of F^, Cl^^, and Fe^^ abundances for K648 are re¬ 
ported here for the first time. The ionic abundances were 
calculated by solving the statistical equilibrium equations for 
more than five levels with the relevant Te and n^, except for 
Ne"^, where we calculated the abundance using a two-energy 
level model. The Fe“^ abundances were solved using a 33- 
level model (from to ^^ 2 )- For each ion, we used the elec¬ 
tron temperatures and densities determined using CEL plasma 
diagnostics. The adopted and for each ion are listed in 
Tabled 

The ionic abundances are listed in Table [T] The final col¬ 
umn shows the resulting ionic abundances, together 

with the relevant errors, including errors from line-intensities, 
electron temperature, and electron density. The ionic abun¬ 
dance and the error are listed in the final row for each ion. 
These data were calculated based on the weighted mean of 
the relevant line-intensity. 

In calculation of the C'*' abundance, we subtracted con¬ 
tamination from [O III] A 2321A to [C ll]A2323 A based 
on the theoretical intensity ratio [Olll] /(A 2326)//(A4363) = 
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0.236. As described above, we did not remove the respec¬ 
tive contributions from and to the [Nil] A5755 A 
and the [O III] A 4363 A line intensities. To determine the 
final abundance, we excluded data determined using the 
[On] AA7320/30A lines. 

We determined the Ne"^ abundance of 2.01(-5) using 
the [Ne ll] A12.80 pm line, which is slightly larger than 
iBover etd] (120061 1.53(-5)). This small disagre ement is ex¬ 
pected to be mainly due to the adopted H/3 flux. iBover et^ 
(l2006h calculated the H/3 flux using the measured HI lines 
at 7.47 pm and 12.37 pm, using the theoretical ratios of 
Hi /(A7.47pm, 12.37pm)//(H/3) with Case B. Their result¬ 
ing /(H/3) was 1.52x10“'^ erg s“^ cm“^. While, we used 
the HST/F656N band-pass flux intensity and corresponding 
HDS spectral scan to sc ale the inten s ities, and And /(H/3) = 
1.07x 10“'^ erg s“' cm~^. IBover et al.l (1200^ used 7): = 10000 
K and = 1700 cm“^ for the Ne'*' and calculations. Our 

plasma diagnostics showed that 10 000 K is low for where 

we used 12 350 K. 

The S^'*' abundance of 2.17(-7) determined using the two 
MIR [S III] lines is approximately the same as that cal- 
culat ed from rSll l] A631 2 A, and is in good agreement 
with iKwitter et al.l (l2003h . who calculated 1.99(-7) using 
/([S III] A 9532 A) = 3.8. However, there was poor agreement 
in the S^'*' abun d ance b etween the most recent m easurements 
by IBover et al.l (120061 2.55(-8)) and our data. IBover et pi 
(120061) calculated the S^'*' abundance using /([S III] A 9532 A) 
= 0.76 measured by iBarkerl (119831) . because they used the 
Spitzer SL module spectra only, where no MIR [S III] lines 
appear. We can exclude the possibility that the discrepancy 
in the S^'*' abundance is due to the flux measurements of our 
MIR [S III] and the choice of /]:. If our flux measurements of 
the MIR [S III], [S III] A 6312 A and the H/3 lines and the 7/ 
selection were incorrect, the abundances from two MIR 
[S III] lines would not match that from [S III] A 6312 A. The 
fine-structure lines are much less sensitive to the electron tem¬ 
perature compared with the other transition lines. The auroral 
lines, such as [S III] A 6312 A, were dependent on the elec¬ 
tron temperature (i.e., the S^'*' abundance determined from 
[S III] A 6312 A is largely dependent on 7/). Our calculated 
abundances from these three were consistent with each 
other, indicating that our flux measurements of the MIR [S III] 
and H/3 lines and the choice of for the S^'*' (and possibly 
also Ne"^ and S^'*') were appropriate. Theref ore, the large dis¬ 
crepancy in S^'*' between bover et ^ (|2006^ and our data may 
have been due to the [S III] A 9532 A flux that was used. 

It is interesting to note the detection of single iso¬ 
tope line candidates [F ll] AA4789.45/4868.99 A, as 
shown in Fig. [T] Together with and ^^Ne, '®F is 
synthesized in the He-rich intershell during the TP-AGB 
phase, and is an «-capture element. The observed [Fll] 
7(A4789.45)/7(A4868.99) of 3.72±0.44 is in agreement with 
the theoretical value of 3.20 calculated using = 10 380 
K and = 3430 cm“^. We excluded the other candidate 
CivA 4789.65A because no Civ lines were detected (e.g., 
CIV A 5801.35 A). Therefore, we conclude that the lines at 
4790 and4869 A are the [Fll] AA 4789.45/4868.99 A, respec¬ 
tively. T he detection of F lines is very rare in Galacti c 
PNe re.g.. lOtsuka eTIfllTtMIZhang & Tiull2?MirhII[T^ . 
Among halo PNe, K648 is the third c ase of suc h F line 
detection rep orted to date; NGC4361 (iLiul 119981) . BoBnl 
(lOtsuka et al.l l2008h . and K648 (this work). We discuss 
whether these lines are [Fll] AA 4789.45/4868.99 A using a 
theoretical model later in the paper. If the two lines do not 


Table 8 

Ionic abundances from RLs. 


x™+ 

Alab 

Multi. 

/(Alab) 

Xm+/H+ 

He+ 

5875.62 A 

Vll 

1.48(+1)±1.38(-1) 

1.02(-l)±6.69(-3) 


4471.47 A 

V14 

4.91(0)±2.78L2) 

9.86(-2)±6.05(-3) 


6678.15 A 

V46 

4.11(0)±5.48(-2) 

9.90(-2)±6.48(-3) 


4921.93 A 

V48 

1.29(0)±5.28(-3) 

9.54(-2)±5.92(-3) 


4387.93 A 

V51 

5.17(-l)±1.04(-2) 

8.34(-2)±6.66(-3) 

1.00(-l)±6.49(-3) 

C2+ 

6578.05 A 

V2 

6.92(-l)±1.08(-2) 

8.42M)±1.33M) 


4267.18 A 

V6 

7.26(-l)±1.33(-2) 

7.32(-^)±9.10C5) 


6151.27A 

V16.04 

4.50(-2)±2.89(-3) 

1.04(-3)±1.31M) 


6462.04 A 

V17.04 

9.81(-2)±8.62(-3) 

9.67(-^)±1.63M) 

8.04M)±1.15M) 


6727.48 A 

V3 

3.71(-2)±2.74(-3) 

2.05M)±1.48(-5) 


6742.15 A 

V3 

4.14(-2)±3.78(-3) 

2.75M)±2.52(-5) 


6744.39 A 

V3 

6.05(-2)±2.80(-3) 

2.87(-^)±1.37(-5) 

2.62(-4)±1.74(-5) 

0^+ 

4641.81 A 

VI 

3.34{-2)±3.45(-3) 

1.21M)±1.64(-5) 


originate from the F^ ion, the prediction cannot fit the fluxes 
of the two lines simultaneously. 


3.1.4. RL ionic abundances 

The RL ionic abundances are listed in Table 0 As we de¬ 
tected C II,III and On RLs, we can compare the elemental 
C and O abundances determined using RLs with those from 
CELs in K648. 

In the abundance calculations, we used the Case B assump¬ 
tion for lines with levels that have the same spin as the ground 
state, and the Case A assumption for lines of other multiplici¬ 
ties. In the final line of each ion series, we give the ionic abun¬ 
dance and the error estimated using the line intensity weighted 
mean. As the RL ionic abundances were not sensitive to the 
electron density with <10* cm“*, we used the atomic data in 
the case of = 10"^ cm“* for all lines. To calculate the He"^ 
abundances, we used r£(He l) = 6710±350K, and the average 
of all Te/He l) data listed in Table|5] except for r^/He l), where 
we used the He 17(A5876)/7(A4471) ratio, which was smaller 
than the other data. We used the Te/BJ) to calculate the 
and abundances. 

We used the multiplet VI Oil A4641.81 A line only, be¬ 
cause the observed HDS spectra were partially contami- 
nated by the absorpt ion lines of the CSPN. According to 
iPeimbert et alJ (120051) . the upper levels of the transitions in 
the V1 OII line are not in local thermal equilibrium (LTE) 
for «£ <10000 cm“*. As the value of calculated using 
the Balmer decrement method was 7500-10000 cm“*, we 
applied the non-LTE corrections using Eq^uations (8)-(10) in 
IPeimbert et al.l (120051) with = 7500 cm“ . 

3.1.5. Nebular ICF abundances 

To estimate the elemental abundances in the nebula, it is 
necessary to correct the ionic abundances that are unseen be¬ 
cause of their faintness or because they lie outside the data 
coverage. We used an ionization correction factor, ICP(X), 
which was based on the IP. The ICP(X) for each element is 
listed in Table |9] The ICP(X)s based on IP are known to be 
inaccurate, particularly in some cases such as N. 

The elemental abu ndances of the nebul a are listed in Ta- 
blefTOl W e referred to lAsplund et al.l (l2009h for N and Cl, and 
iLoddersI (120031) for the other elements. 

The RL C abundance was almost identical to that of the 
CEL C, that is, the C abundance discrepancy factor (ADE) 
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Table 9 

The ionization correction factors (ICFs). 


X 

Line 

ICF(X) 

X/H 

He 

RL 

S++S2+ 

ICE(He)He+ 

C 

CEL 

{■S').. 

C++ICE{C)C2+ 


RL 


ICE(C)C2++c3+ 



V C-+ /cel 

N 

CEL 

(^)cEL 

ICE(N)N+ 

O 

CEL 

1 

0++0-+ 


RL 

(i^)cEL 

ICE(0)02+ 

F 

CEL 

(^)cEL 

ICE(F)E+ 

Ne 

CEL 

1 

Ne++Ne2+ 

S 

CEL 

1 

S+ + S2+ + S^+ 

Cl 

CEL 

(^) 

ICF(C1)CP+ 

Ar 

CEL 

s 

w 

ICF(Ar)Ar2+ 

Ee 

CEL 

(^)cEL 

ICF(Fe)Fe2+ 


n(C)RL/n(C)cEL = 1.17±0.75, whereas the O ADF was large, 
«(0)rl/«(0)cel = 2.99±0.55. The RL C abundance is greater 
than the RL O abundance. The RL C/O ratio of 17.46±7.07 
agrees with the CEL C/O ratio of 6.83±3.63 within error. The 
(C/0)rl/(C/0)cel ratio is 2.56±1.71. It follows that these 
C/O ratios indicate that K648 is a C-rich PN. 

The aforementioned O ADF value in K648 is approxi¬ 
mately the same as the O^"^ ADF = 2.99±0.46. The O^"^ 
ADF has be en reported for other C- rich halo PNe, i.e., H4-1 
(1.75±0.36. | Otsuka & TaiitsuHMTl and BoBnl (3.05±0.54, 
lOtsuka et^ 120101) . The smaller ADF of the O^"^ found in 
H4-1 ma y be d ue to temperature fluctuations proposed by 
iPeimberj (119671) : however, the relatively large ADF of O^"^ 
found in K648 is too large to be explained by temperature 
fluctuations. Therefore, as with BoBnl, we should seek other 
plausible solutions to ex plain the O^'*' a bunda nce discrepancy 
in K648. For BoBnl. TOtsuka et al.l (1201 Oh suggested that 
the bi-abundance pattern may solve the and Ne^'*' abun¬ 
dance discrepancy. The RL abundances in the nebula may 
correspond to abu ndances in the stellar wind, as seen in the 
C-rich PN IC418 (iMorisset & Georgievll2009l) . It should be 
no ted that the stellar C and O abundances for K648 examined 
by iRauch et'dl (l2002h are closer to our RL C and O abun¬ 
dances (C = 9.00 and O = 9.00). In the following section, we 
determine the stellar abundances of K648 using the FUSE, 
HST/COS, and HDS spectra, and check for correlations with 
the stellar abundances of the nebular RL C and O values in 
K648. 

Table [TT] lists the nebular elemental abundances of K648. 
These data were det ermined using the se mi-e mpirical ICF 
metho d, except for [Howard et al.l (119971) and lAldrovaridil 
(Il980h . who obtained the abundances using photo-ionization 
(P-I) models. We determined the abundances of Ne, S, and 
Ar, as well as that of CEL C, and added those of RL O, and 
the CEL F, Cl, and Fe using the HDS and Spitzer/IRS spec¬ 
tra for many ionization stages. Our measurements show good 
agreement with those reported previously, with the exception 
of those for C and N. Scatter in the CEL C abundance may 
be due to the use of Tg for the C^'*' abundance and/or the H/) 
flux measurements, because the emissivity of the C III] lines is 


very sensitive to Tg. Note that the observation window of the 
international ultraviolet explore (lUE) is very large for K648 
(window dimension: 10.3x23 arcsec^ elliptical shape). The 
scatter of N abundance may be due to the use of ICF(N). We 
will ch eck t he CEL C and N abundances in the P-I model in 
Section iTAl The F abundance is comparable to that in BoBnl 
(F ^ = 5.98.IC>tsuka et al.ll2010l) . The Ne abundance reported 
bv iBover et al.l (120061) was performed by adding the Ne'*' abun¬ 
dance determined from the Spitzer/IRS spectrum, whereas 
others did not calculate the Ne'*' abundance. For this reason, 
our Ne ab undance is larger th an has been reported previously, 
except for lBover et al.l (|2003). 

3 . 2 . Absorption line analysis 

We employed a spectral synthesis fitting method to inves¬ 
tigate the elemental abundances in the photosphere of the 
CSPN of K648 using O- type star grid models (OStar2002 
grid) based on TLUSTY (iLanz & Hubenvl2003l) . which con¬ 
siders 690 metal line-blanketed, non-LTE, plane-parallel, and 
hydrostatic model atmospheres. We considered the 8 ele¬ 
ments He, C, N, O, Ne, S, P, S, Fe, and Ni, together with 
approximately 100000 individual atomic levels from 45 ions 
(see Table 2 of iBouret et al.ll200^ . 

3 . 2 . 1 . Modeling process 

We found [Ar,Fe/H] abundances of -1.96 and -2.45 from 
the nebular line analysis, respectively, and sorted models with 
a metallicity of Z = 0.01 and 0.001 Zq from the OStar2002 
grid models. All of the initial abundances in these models 
(except He) were set to [X/H] = -2 (0.01 Zq) models and -3 
(0.001 Zq). The initial ratio of He/H abundances was set to 
0.1 in both models. _ 

Following iBouret et alJ (120031) and iRauch et al.l (120021) . we 
determined logg, Teff, and the He/H abundance ratio, which 
are the basic parameters used for characterizing the photo¬ 
sphere. First, we generated models with [X/H] = -2.3, cor¬ 
responding to 0.005 Zo, by interpolating between the 0.01 
and 0.001 Zq grid models using the IDL programs INTRP- 
MOD and INTRPMET. We set the microturbulent velocity 
to 5 km s“' and the rotational velocity to 20 km s“*, because 
models with these values were found to fit the absorption line 
profiles in the FUSE and the HST/COS spectra, as well as 
the HDS spectrum. Before attempting to determine Teff and 
log g using the stellar absorption lines, we ran the SED mod¬ 
els using CLOUDY (iFerland et al.lll998h to find the ranges 
of Teff and logg in the TLUSTY models. These CLOUDY 
and SED models maintain the initial photosphere abundances 
(i.e., He/H = 0.1 and 0.005 Zq). We found that the models 
can reproduce the observed HST/WFPC2 F547M flux density 
and the emission line fluxes if we use the incident SED gener¬ 
ated using the TLUSTY model with 0.005 Zq, 7’eff'^34000- 
40000 K, and logg~3.5-4.1 cm s“^. 

Using the 0.005-Zq grid models, we determined logg 
and He/H by monitoring the chi-squared value of the HDS 
He IIA 4541 A and the synthesized line profiles of this line. 
We ran grid models with T^ff = 35 000-41000 K (in 100-K 
steps), logg= 3.5-4.1 cm s“^ (inO.Ol-cm s“^ steps), and He/H 
= 10.98-11.06 (in steps of 0.01). We used SYNSPEC to gen¬ 
erate synthesized spectra. We set the spectral resolution toR = 
33 500 and used a heliocentric radial velocity of -125.30 km 
s“^ determined using the He IIA 4541 A absorption line before 
running SYNSPEC. We monitored the spectrum in the range 
4535-4547 A. The best fit was given by logg = 3.96±0.02 cm 
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Table 10 

The elemental abundances from CEL and RLs. 


X 

Types of 
Emissions 

X/H 

log(X/H)+12 

[X4I] 

log(Xo/H)+12 

ICF(X) 

He 

RL 

L04(-l)±6.82(-3) 

11.02±0.03 

+0.09±0.03 

10.93±0.01 

1.04±0.01 

C 

CEL 

9.41(-4)±3.75M) 

8.97±0.17 

+0.58±0.18 

8.39±0.04 

1.33±0.24 

C 

RL 

L10(-3)±5.54M) 

9.04±0.22 

+0.65±0.22 

8.39±0.04 

1.04±0.67 

N 

CEL 

2.28(-6)±5.35(-7) 

6.36±0.10 

-L47±0.11 

7.83±0.05 

4.02±0.81 

O 

CEL 

5.39(-5)±3.84(-6) 

7.73±0.03 

-0.96±0.06 

8.69±0.05 

1.00 

O 

RL 

L61M)±2.72(-5) 

8.21±0.07 

-0.48±0.09 

8.69±0.05 

1.33±0.13 

F 

CEL 

2.60(-7)±6.70(-8) 

5.42±0.11 

+0.96±0.13 

4.46±0.06 

4.02±0.8I 

Ne 

CEL 

2.75(-5)±2.19U6) 

7.44±0.03 

-0.43±0.11 

7.87±0.10 

1.00 

S 

CEL 

2.53(-7)±3.93(-8) 

5.40±0.07 

-L79±0.08 

7.19±0.04 

1.00 

Cl 

CEL 

3.76(-9)±1.28(-9) 

3.58±0.15 

-L92±0.33 

5.50±0.30 

1.19±0.29 

Ar 

CEL 

4.00(-8)±1.17(-8) 

4.60±0.13 

-L95±0.15 

6.55±0.08 

1.19±0.29 

Fe 

CEL 

L06(-7)±2.84(-8) 

5.02±0.12 

-2.45±0.12 

7.47±0.03 

4.02±0.81 

Note. 

— The types of emission 

line used to calculate the abundances 

are shown in the second column, the number densities of each element relative to hydrogen are 

listed in the third 

column, the fourth column lists the number densities, where log^Q /i(E) 
list the solar abundances and the ICF values that were used. 

= 12, the fifth column lists the logarithmic number densities relative to the solar value, and the final two columns 


Table 11 

Comparison of nebular elemental abundances. 

References 

He 

C 

N 

0 

F 

Ne 

S 

Cl 

Ar 

Fe 

This work (RL) 

11.02 

9.04 


8.21 







This work (CEL) 


8.97 

6.36 

7.73 

5.42 

7.44 

5.40 

3.58 

4.60 

5.02 

Bover et^^(2W^ 






738 

T53 




Kwitter ^^^^^31 

11.00 


6.48 

7.85 


7.00 

5.30 


4.60 


Howard et al. (1997)“ 

Henrv et^^_^96)'’ 

10.98 

8.50 

6.72 

7.61 


6.57 

6.11 


3.72 


10.92 

8.29 

6.66 

7.62 


6.47 





Adams et al. Q984I'’ 

11.02 

8.73 

6.50 

7.67 


6.70 





A Idrovan^^^^OTr 

10.90 

8.45 

6.37 

7.53 


6.40 

5.60 




Torres-Peimbert & Peimbert (1979) 

10.99 


<6.39 

7.82 


6.79 

<6.22 


<5.52 


Hawlev & Miller (1978) 

11.00 


7.11 

7.65 


6.40 






^ from the photo-ionization models. 

^ The C abundance is from CEL C lines. 


and He/H = 11.05±0.02. In this process, we estimated 
Tpff =37 000 K. Ou r data are in good agreement with those of 
iRauch et~an (I2002h . who reported logg = 3.9±0.3 cm s“^, 

= 39000±2000 K and He/H = 10.9±0.3 obtained using their 
non-LTE model. 

We determined Eeff and the abundance of C assuming that 
logg = 3.96 cm s“^ and He/H = 11.05. Here, we used 
the C III and C IV lines in the FUSE and the HSTJCOS 
spectra, including C III A 1246/47 A, C IV A 1107/08 A and 
CIV A 1230/31 A. At approximately Eeff = 35 000-41000 K, 
the strengths of the C III lines were sensitive to Eeff, whereas 
those of the C IV lines were not. Therefore, we can determine 
Teff accurately and the abundance of C simultaneously using 
a plot of the C abundance as a function of Teff. We find Teff = 
36 360±700 K and C = 9.38±0.10. 

Using logg = 3.96 cm s“^ and Teff = 36 360 K, we de¬ 
termined the N, O, Ne, P, and Fe abundances to match the 
observed line profiles. We used SPTOOPQfor line identifi¬ 
cation. The N abundance was obtained using N III A 1243 A 
and N IV A 1719 A. The O abundance was found from the 
many O III lines around 3774 A in the HDS spectrum and at 
AA 1149/51 A, O IV AA 1342/44, and O v AA 1371 A. The Ne 
abundance was found from Ne III A 1257 A only. The P abun- 

^ SPTOOL is a software package for analyzing high-dispersion stel¬ 
lar spectra (i.e., line identification, determination of radial velocity, inves¬ 
tigation of the atmospheric parameters, such as turbulent velocities or ele¬ 
mental abundances), developed by Youichi Takeda. We also used the AT- 
LAS9/WIDTH9packages written by R. L. Kurucz. 


Table 12 

The central star properties determined using the TLUSTY model. 


Parameters 

Values 

Basic Parameters 

Tcff 

36 360±700 K 

log? 

3.96±0.02 cm s^^ 

Photosphere Abundances (log(H) = 12) 

He/H ([He/HJ) 

1L05±0.02 (+0.1240.02) 

C/H ([C/H]) 

9.3840.02 (+0.9940.04) 

N/H ([N/H]) 

6.5340.10 (-1.3040.11) 

0/H ([O/H]) 

8.3640.10 (-0.3340.11) 

Ne/H ([Ne/H]) 

8.2140.10 (+0.3440.14) 

P/H ([P/H]) 

3.6440.10 (-1.8240.11) 

Fe/H ([Fe/H]) 

5.2340.10 (-2.2440.10) 


dance was determined from the P y AA 1118/28 A, and the Fe 
abundance from Fe V AA 1448/56 A. 


3.2.2. Comparisons between stellar and nebular abundances 

The resulting spectrum synthesized using the TFUSTY 
and the observed FUSE and HST/COS spectra are shown in 
Fig. 0 The parameters, including the elemental abundances, 
are listed in Table [T2l The derived stellar abundances are also 
listed in Table fTTI for comparison. We were unable to detect 
any F absorption lines, e.g., F V A 1082/87/88 A due to the low 
S/N ratio. The detection of the single isotope ^'P is interest¬ 
ing because phosphorus (along with fluorine) is an «-capture 
element that is synthesized in the He-rich intershell during the 
TP-AGB phase. 
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Figure 8. The line-profiles of the selected lines observed in the FUSE and 
COS spectra (black lines) and the synthesized spectrum calculated using the 
TLUSTY model (red lines). The wavelength is shifted to the wavelength in 
vacuum. 


We found that the stellar C and O abundances were close 
to the nebular RL abundances; however, the stellar He and 
Ne abundances were larger than the nebular abundances. The 
stellar N and Fe abundances were comparable to the nebular 
abundances. We may expect slightly higher stellar C, O, and 
Ne abundances than the nebular abundances, as the former are 
indicative of more recent products of AGB nucleosynthesis. 
These three elements are synthesized in the He-rich intershell 
during the AGB phase, and are then brought up to the stellar 
surface via the TDU. Note that the stellar C/O and the Ne/O 
ratios (10.75±2.43 and 0.75±0.23, respectively) are in good 
agreement with nebular ratios C/O (17.46±7.07 in RL and 
6.83±3.63 in CEL) and Ne/O (0.51±0.06) in CEL. Although 
it is difficult to determine whether the RL or CEL abundance 
represents the nebular C and O chemical abundances in K648, 
the similarity of the C/O ratios determined from the RLs and 
CELs indicates a positive correlation with the stellar abun¬ 
dance. 

K648 shows large stellar and nebular CEL [O/Ee] abun¬ 
dances (L91±0.15 dex versus L49±0.13 dex). The [Ne/Ee] 
abundances were also large (2.01 ±0.16 dex in the CEL and 
2.58±0.18 dex in the stellar region). It has been reported 
that metal-poor stars in the Milky Way exhibit large [a/Ee] 
abundances, where the a-elements include O, Ne, Mg, Si, and 
Ca. The effect is greatest for the most metal-poor populations, 
such as members of the stellar halo and, in particular, in the 
rO/Ee l (see, for example. lMcWilliamlfTo^lFeltzing & Chibal 
l2013h . This is interpreted as a consequence of time delay in Fe 
production from Type la SNe relative to the a-elements from 
core-collapse SNe. The a-elements are mainly produced by 
Type II SNe. Both types of SNe should produce Fe in the 
proportio ns of ^^1/3 for Type II and ~2/3 for Type la SNe. 
For M15, ISobeck et al.l (1201 II) reported that three red giant 
branch (RGB) stars (([Fe/H]) = -2.55), exhibited O abun¬ 
dance of 6.75-7.03 and the [O/Fe] of h- 0.62 to h- 0.85 (([O/Fe]) 
= h- 0.75). If we observe an RGB star with [Fe/H] = -2.3, 
[O/Fe] should be -1-0.50, which corresponds to an O abun¬ 
dance of 6.89. 


The difference between the [O/Fe] of M15 RGB stars re¬ 
ported bv ISobeck et al.l (1201 Ih and that of K648 suggests that, 
in K648, O synthesis was >0.9 dex during the TP-AGB phase. 
The TP-AGB phase nucleosynthesis process can contribute 
to enhancement of O and Ne abundances in the helium con¬ 
vective zone with formed from mixed protons as an n- 
source using a nuclear network from H through S. The abun¬ 
dance of '®0 may increase in proportion to the square root 
of the amount of mixed until it reaches a significant frac¬ 
tion of '^C, whereas the abundance of ^^Ne may increase in 
proportion to t he amount of mixed and a ttains half of 
the mi xed (iNishimura et alJl2009t) . Indeed. ILugaro et aP 
(1201 2h demonstrated that [Fe/H] = -2.19 AGB stars c an sy n- 
thesize significant quantities of O and Ne (see Section l4Tt . 

3 . 2 . 3 . The core-mass of the CSPN 

The core-mass of the CSPN can impose a significant con¬ 
straint on the initial mass of the progenitor. Through construc¬ 
tion of a TLUSTY model atmosphere, we obtained the H\ 
spectrum of the stellar photosphere. Using the Hx and the ob¬ 
served //5r/WFPC2 F547M flux density Ix listed in Tabled] 
we de termined t he cor e-mass of the CSPN Me using Equation 
(1) of lShipmanl(ll979l) . i.e., 

Ix=ATtHxR^D-^, (3) 

g = GMeR-\ (4) 

where R is the radius of the CSPN, D is the distance to K648 
from us, g is the surface gravity of the CSPN, and G is the 
gravitational constant. 

Using the synthesized spectrum from TLUSTY model at¬ 
mosphere fitting, we found that Hx = 4.53(-i-7) erg s“' cm““ 
A“' at A 5483.88 A by taking the transmission curve of the 
HST WFPC2/F547M band into account. Recent measure- 
met^ of the distance t o M15 have been repor ted by iReidI 
(fT^ 1 2.3±0.6 kpcL iMcNamara etalJ (l2(M 9.98±0.47 
kpc), and ivan den Bosch et al.l (120061 10.3±0.4 kpcl. We find 
logg = 3.96±0.02 cm s“^, as determined in Section U.2.11 

If we use the average distance amongst these distance mea¬ 
surements, i.e., 10.9±0.5 kpc. we obtain = 0.68±0.07 Mq 
andR= L43±O.O8R0 using Equations (O and ([Hi. Using the 
most recent data, i.e., D = 10.3±0.4 kpc, the values of Me and 
R are Me = 0.6 1 ±0.06 an d L35±0. 08 Rrx, whic h are in 
agreement with iBianchi et al.l (1200 Ih and iRauch et al.l ([2002) ■ 
We found that R = 1.3 Rq and Me of O.62±O.1O_M0 with D 
= 10.3 kpc and logg = 4.0 cm s“^. [Ranch et aP (l2002h calcu¬ 
lated Me = 0.57 Mq from the theoretical Teff-logg diagram. 
The exact value of the Me is still dependent on the choice of 
dist ance. We will discuss the initial mass of K648 in section 

lO 


3.3. Dust features in the Spitzer//R5 spectrum 

As discussed in Section 12.71 K648 exhibits the 6-9 p,m 
PAH band, the 11.3 pm PAH band, and the broad 11 pm fea¬ 
ture. These PAH bands are some times seen in C-rich PNe , 
such as BD-(-30° 36 39 (C/O = L59. lBernard-Salas et al.ll200l 
[Waters et al.lll998h. as well as O-rich PNe such as NGC6302 
(C/O = 0.43. [Wright etaPIMll: [Molster et akllMll . Both 
BD-i-30° 3639 and NGC6302 exhibit strong crystalline sili¬ 
cate features at 23.5, 27.5, and 33.8 pm, which have never 
been observed in K648. In addition, the 9 and 18 pm fea¬ 
tures attributed to the amorphous silicate were also not seen 
in K648. Therefore, we concluded that K648 is a C-rich gas- 
and-dust PN. 




































































12 


Otsuka et al. 



Figu re 9. The 5-15 /im spec tra of C-rich halo PNe K648 (black), H4-1 (blue, 
from ITaiitsu & Otsukal2Q14D . andBoBnl (red, from lOtsuka et alJ2QlCD . The 
local continuum was subtracted from the observed flux, which was then nor¬ 
malized to that at 8.6 /im. The data for K648 are obtained from from the SL 
module. The 5-10/im spectra of H4-1 andBoBnl were obtained from the SL 
module and the remaining data from the SH module. 


Table 13 

The results of Gaussian fittings to the PAH bands, the possible 6.85 ^m 
aliphatic feature, and the broad 11 pm band. 


Ac 

(/xm) 

FWHM 

(/xm) 

F{X) 

(erg cm“^) 

/(A) 

[/(H/3) = 100] 

6.25±0.01 

0.22±0.02 

2.25(-14)±1.88(-15) 

2.20±0.21 

6.47±0.0P 

0.13±0.03 

4.24(-15)±9.82(-16) 

0.41±0.10 

6.85±0.02'’ 

0.21±0.04 

6.01(-15)±1.65(-15) 

0.59±0.16 

7.83±0.01 

0.25±0.03 

1.34(-14)±2.12(-15) 

1.31 ±0.22 
2.17±0.13‘‘ 

8.73±0.01 

0.16±0.04 

5.37(-15)±1.55(-15) 

0.53±0.15 

11.31±0.01 

0.25±0.01 

1.00(-14)±6.46(-16) 

0.98±0.08 

11.81±0.03" 

1.97±0.14 

1.23(-13)±1.09(-14) 

12.10±1.21 


^ He I 6.47 pm. 

^ The complex o f the p ossible aliphatic 6.85 ^m and the [Aril] A 6.99/im 
lines. See Section [3.3.2l regarding the respective intensities. 

^ The broad 11 pm band. 

^ The extrapolated value using the PAH /(A 7.7 pm)/l{\^.6 pm) ratio in 
BoBnl. 

Figure |9] shows the 5-15/xm spectrum, where the local 
dust continuum was subtracted by fourth-order spline fit- 
tin g, using the s ame t echnique as applied for C-rich PNe 
by lOtsuka et alJ (12014) . The flux density was then normal¬ 
ized to the intensity of the 8.6 fim PAH band. For compari¬ 
son, we als o show the Spitzer/IRS spectra of the C-rich halo 
PNe H4-1 (iTaiitsu & Otsuk^lmj) . as well as that of BoBnl 
(lOtsuka et al.ll2010l) . We discuss the dust features in more de¬ 
tail below. 

3.3.1. The 6-9 fim and 11.3 fxin PAH bands 

The 6-9 /xm and 11.3 /xm PAH band profiles are remarkably 
similar to those of H4-1 and BoBnl, although the intensity 
peak of the 7.7 pm PAH in K648 is smaller than those of 
BoBnl and H4-1, which is attributed to noise around 7.7 pm. 

We measured the central wavelength A^, FWHM, flux F(X), 
and relative intensity /(A) of each PAH band by single Gaus¬ 
sian fitting, and the results are shown in Table [13] For the 
6.2 pm band, we employed a double Gaussian component fit, 
where one component corresponds to the 6.25 pm PAH band 
a nd the other to the He IA 6.47 pm. 

iPeeters et al.l (l2002h examined the profiles of the 6.2, 7.7, 
and 8.6 pm PAH bands using ISO/S'WS spectra, and classified 


the spectra into Classes A, B, and C according to the peak po¬ 
sitions of each PAH feature. Class B PAHs are frequently seen 
in C-rich PNe, including BoBnl and H4-1, and have a peak in 
the range 6.235-6.28 pm, a stronger component at ^7.8 pm 
than at 7.6 pm, and a peak at >8.62 pm. The 6.2, 7.7, and 
8.6 pm features in K648 satisfy the definition of a Class B 
PAH spectrum. 

According to the classification of the 11.3 pm PAH profiles 
by Ivan Diedenhoven et al.l (12004) . the 11.3 pm PAH in K648 
falls under Class Bn 2 , with a peak at ^11.25 pm. Many C- 
rich PNe in the Magellanic Clouds also have this class of I’AH 
(iBernard-Salas etal1l2009^ . 

3.3.2. The 6.85 fim aliphatic feature? 

K648 exhibits a weak broad feature at 6.85 pm, which 
might be a combination of the 6.85 pm aliphatic fea- 
ture (CH 2 r asymm etric deformation) and [Aril] A6.99pm. 
lOtsuka et al.l (12014) established a relationship between Teff 
and the /([Arm] A8.99pm)//([Arll] A6.99pm) ratio in C- 
rich PNe based on P-I models with CLOUDY co de. U sing their 
Equation (Al) and Teff = 37 100 K (see Section iTTl l. we found 
that the 6.85 pm aliphatic feature and the [Ar ll] A 6.99 pm in¬ 
tensities are 0.37±0.09 and 0.21±0.13, respectively, where 
the H/3 inten sity is 100. _ 

Following iLi & Drainel(l2012^ . we estimated the number ra¬ 
tio of C-atoms in aliphatic form relative to those in aromatic 
form using the 6-9 pm PAH band, i.e., Ac,aiiph/Ac,arom- As 
we underestimated the 7.7 pm PAH flux in K648, we extrap¬ 
olated a 7.7 pm PAH intensity of 2.17±0.13 using the PAH 
7(7.7pm)//(8.6pm) ratio of 4.11±0.13 measured in BoBnl. 
Our derivation is Ac,aiip[i/Ac,arom of ~0.1-0.4, indicating that 
<29% of the C-atoms exists in aliphatic form in K648. 

For a more accurate estimate of the number of C-atoms in 
the aliphatic form, L-band spectroscopy is useful to check 
for the existence of the 3.4 pm aliphatic feature, as well as 
7(3.3 pm PAH)/7(3.4pm aliphatic). 

3.3.3. The broad 11 pm band 

K648 exhibits the broad 11 pm feature, which is fre¬ 
quently seen_jn_Galactic an d Magellanic C-ri ch PNe 
le.g.. lO tsuka et "^ 12014 iBernard-Salas et aTI 120091 
IStanghellini et al.r^lilL The band profile appears to 
show an almost flat portion in the range 11.4-12.2pm. 
However, as shown in Fig. [TOl al. the resulting band profile 
did not exhibit a flat top after removal of the 11.3 pm PAH 
band and the atomic lines. 

The results of Gaussian fitting are also listed in Table [13] 
The FWHM of the 11 pm band is comparable to those for 
H4-1 (2.08±0.05 pm) and BoBnl (1.85±0.39pm); however, 
Ac was slightly blue-shifted (12.28 ±0.06 pm in H4-1 and 
12.30±0.08 txm in BoBnl ). Our results corroborate those of 
IBernard-Salas et aTI (120091) . who reported that the profile and 
the central wavelength of the 11 pm band in MC PNe differ 
from source to source. 

There is some debate regarding the origin of the broad 
11 pm feature. Silicon carbide (SiC) is one possible ex- 
planation for the feature at 11 pm in C-rich MC PNe 
(IBernard-Salas et aDl2009^ . In Fig. [T0[al. as a SiC template, 
we show a comparison with the 11 pm band profile of the 
Galactic solar metallicity C-rich AGB star W Ori, extracted 
from the arc hive ISO/S'WS spe c trum. These data w ere down¬ 
loaded from ISloan et al.l (120031) . lAbia et all (120021) reported a 
C/O ratio of 1.005, and a metallicity of [M/H] = -(-0.05. The Ac 
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Figure 10. The broad 11 /im band profiles. The 11.3 PAH and the atomic 
lines were subtracted out via Gaussian fitting in both (a) and (b). The flux 
density is shown nonnalized to the intensity peak. The spectral resolution 
of the C-rich AGB star W Ori and the proto PN IRAS Z02229+6208 was 
adjusted to match that of K648. 

(11.2 /j,m) and FWHM (1.51 /J,m) in the 11 //m band of W Ori 
differ significantly from those measured for K648. The 11 fim 
band profile in W Ori may be fitted to an absorption efficiency 
Q\ of a spherical g- SiC grain (or 6- H SiC, hexagonal unit 
cell) calculated from iPegouriel (Il988h . which peaks sharply 
at ^11.2 fim and has an FWHM of ^1.2/xm. Howe v er, we 
must be careful with Ac in W Ori. iLeisenring et alJ (l2008h 
demonstrated how the C 2 H 2 absorption band around 13.7 p,m, 
as well as the SiC self-absorption band around 10 pm affect 
the central wavelength of SiC in AGB stars such as W Ori. 
They argued that the C 2 H 2 absorption band suppressed the 
long-wavelength part of the feature at 11 pm, and caused the 
central wavelength to be blue-shifted. This may be the case 
for W Ori. We could fit neither Ac nor the FWHM of the 
11 pm band , even with a continuous distribution of elli psoids 
(CDE, e.g.. lBohren & HuffmanllMl IMin et aUlMlll of g- 
SiC; we find Ac ^ 11.8pm and a FWHM of ^2 pm using the 
Q\ for the CDE g-SiC, as shown by the by the gray line in 
F ig-IMa). _^ 

iKwok et al.l (120011) argued that a collection of out-of-plane 
bending modes of aliphatic side groups attached to an aro¬ 
matic ring could result in to the broad 11 pm feature. In¬ 
deed, we found the feature at 6.85 pm in K648, correspond¬ 
ing to possibly aliphatic C. Eigure fTOl' bl shows a compari¬ 
son of the 11 pm band profiles of K648 and the proto-PN 
IRAS Z02229H-6208. The [C/H] and [M/H] ab undances of 
this p roto-PN are -1-0.29 and -0.50, respectively, (iReddv et alJ 
Il999h . The measured values of Ac and the EWHM of IRAS 
Z02229-I-6208 are 11.68±0.02pm and 1.81 ±0.05 pm, re¬ 
spectively. The 11 pm band profile in IRAS Z02229-I-6208 
exhibits a good fit to that of K648, except for th e <11 p m part 
of the 11 pm band. According to iKwok et aH (120011) . K648 



Photon Energy [eV] 

Figure 11. The SED of the CSPN of K648 synthesized using TLUSTY. 
This SED was adopted in the dust+gas SED model using CLOUDY code as 
the incident SED. 



Distance from the CSPN (pc) 

Figure 12. The hydrogen density radial profile used in the CLOUDY P-I SED 
modeling. 

may have a few cyclic alkanes, which contribute to the 9.5- 
11.5 pm part of this band (see Pig. 4 of lKwok et alJLOOHl . 

The low metallicity of K648 implies a very low abundance 
of Si. Indeed, we did not detect any lines corresponding to Si 
in either the nebula or the central star. Therefore, we expect 
that the broad 11 pm band profile in K648 is attributable to a 
wide variety of alkane and alkene groups attached to hydro¬ 
genated aromatic rings, rather than to SiC. 

3.4. Radiative transfer modeling and SED fitting 

We constructed an SED model to investigate the physical 
conditions of the gas and dust grains and derive their masses 
using Cloudy cl0.00. The quantity of dust mass formed in 
extremely metal-poor objects such as K648 is of interest. The 
gas mass as well as the core-mass of the CSPN are required 
to unveil the origin and evolution of K648 via a comparison 
of these parameter values with the results of AGB nucleosyn¬ 
thesis models. 

3.4.1. Modeling approach 

We attempted to fit the observed SED in the range 0.1- 
160 pm, assuming that the dust in K648 is composed of PAH 
molecules and amorphous carbon (AC) grains. No SiC grains 
were considered to fit the broad 11 pm feature. 
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The distance to K648 is required to compare our model val¬ 
ues with the observed fluxes; here we assumed a distance of 
10.9 kpc. For the incident SED from the central star, we used 
the synthesized spectrum of the central st ar of K648 using the 
TLUSTY model, as discussed in Section [T2] The input SED 
is shown in Eig. [TT] We adjusted the input SED to match 
the de-reddened absolute V-band magnitude of -0.528 mea¬ 
sured from HSTfPSAl'bA photometry of the CSPN. The num¬ 
ber of Lyman continuum photons A^Lyc with >13.5 eV was 
6.92 (h-45) s“', as determined from the synthesized spectrum 
of the CSPN. 

The P-I model construction with the CLOUDY and other 
codes involves an ad hoc nebular geometry and central stel¬ 
lar property modification. Until it gives a right prediction to 
the line intensities and continuum, one must adjust not only 
the chemical abundances but also the model nebular geometry 
with a new value close to the observation indication. In order 
to tune the other diagnostically indicated physical properties, 
e.g., electron temperature, one even needs to consider other 
chemical elements which were not observed at all. We em¬ 
ployed the observed values of the gas-phase elemental abun¬ 
dances listed in TablefTOlas initial estimates, and refined these 
to match the observed line intensities of each element. We 
considered the RL and CEL C line fluxes and the observed 
CEL O line fluxes to determine the nebular C and O abun¬ 
dances, respectively. We revised the transition probabilities 
and collisional impacts of C III], [N ll], [O II, III], [Ell,iv], 
[Ne ll,lll,iv], [S ll,lll], [Clll,lll] and [Ar ll,lll,iv], which 
were the same as those used in our semi-empirical ICE abun¬ 
dance calculations (i.e., using IP coincidence method). The 
abundances of other elements were fixed to be constant, with 
[X/H] =-2.3. 

We determined the radial hydrogen density profile of 
the nebula based on the radial intensity profile of the 
//57’/WFPC2 P656N image using Abel transformation, and 
assuming spherical symmetry. We fixed the outer radius to 
Rout = 2.1" (0.11 pc) and the inner radius to Rin = 0.14" 
(0.0072 pc). We used a constant filling factor of e = 0.5. The 
hydrogen density radial profile is shown in Pig. [12] The Rout 
that was used corresponds to the Stromgren radius (0.11 pc), 
assuming = 10"^ K, = 3000 cm“^, and the same 

values of e and ALyc 

We assume that both PAH molecules and AC grains exist in 
the nebula, and that the observed IR-excess from MIR to PIR 
wavelengths is due to the thermal emission from these species. 
We assumed spherical AC grain s and PAH molecule s. The op¬ 
tical con stants were taken fromibraine & lH (l2007h for PAHs 
and from iRouleau & Martini (1 199 ll) for the AC grains. Por the 
PAHs, we assumed that the radius was in the range of 0.0004- 
0.0011 pm (i.e., 30-500 C atoms) with an a~^-^ size distribu¬ 
tion. Por the AC grains, we used th e standard interstella r dust 
grain size distribution reported by iMathis et al.l (|l973), i.e., 
an a~^'^ size distribution, but with a smaller radius of 0.0005- 
0.010 pm, which was determined by running several test mod¬ 
els. 

To evaluate the degree of accuracy of the model fitting, we 
calculated the chi-square (x^) value from the 39 gas emis¬ 
sion fluxes, 10 gas-phase abundances, and the five broad band 
fluxes, as well as the 15 flux densities of the features of inter¬ 
est from UV to FIR wavelengths. 

3.4.2. Modeling results and SED fitting 

Figure [13] shows the predicted SED, the observed spectra, 
and the band flux densities. The predictions were taken at the 


Table 14 

The properties from the CLOUDY P-I model. 


Parameters 

Values 

Central Star 


-U.528, measured from /h54/M obs 

L, 

3076 Lq 

Lff 

36 360 K 

logg 

3.96 cm 

Distance 

10.9 kpc 

Nebula 

Abundances^ 

Heill.OO, C:8.71, N:6.96, 0:7.82, 

(Iogn(X)/H(H)+12) 

F:5.4I, Ne:7.02, S:5.48, Cl:3.44 


Ar:4.44, Fe:5.48, Others:[X/H] = -2.3 

Geometry 

Spherical 

Shell size 

Rin = 0.14 " (0.0072 pc), Rout = 2.1 " (0.125 pc) 

rtH 

See Fig. [12] 

filling factor 

0.50 

Iog/(H/3) 

-11.972 erg s“^ cm“^ (de-redden) 

mg 

4.81(-2)M0 

Dust in Nebula 

Composition 

PAHs, amorphous carbon (AC) 

Grain size 

0.0005-0.010 n-m for AC 


0.0004-0.011/rm for PAH 

U(PAHs) 

140-472 K 

U(AC) 

99-290 K 

mrffTot.)'’ 

4.95(-7) Mq 

mrf(Tot.)/mg 

1.029(-5) 


“ The relative eiTor of the gas-phase elemental abundances is within 0.2 dex. 
The total dust mass of the PAH and AC grains. 


matter-bounded radius near the Stromgren edge (or at the ra¬ 
dius close to the ionization-bounded radius of the P-I model 
nebula). This provides an appropriate level of nebular excita¬ 
tion, e.g., for 0“^/(0'^H-0^'^). Note that the observed and pre¬ 
dicted nebular ra tios O^Y(0'^ + O^'| ' )~0.75 (0.92 in BoBnl 
and 0.67 in H4-L lOtsuka et al1l2010l : IOtsuka & Taiitsull2()T^ 
were large despite the cool CSPN of K648. Such a high ratio 
indicates that K648 could be a matter-bounded nebula, where 
the edge of the mass distribution falls inside the Stromgren 
edge, rather than an ionization-bounded nebula, and also it 
might be related to the small nebula mass. 

The fitted elemental abundances, gas mass nig, dust mass 
nid, and dust temperatures Td are listed in Table [14] The 
third and fourth columns of Table [15] show a comparison 
of the predicted fluxes and flux densities with the observed 
data. The discrepancies of each flux and each flux density 
between the observation and model are listed in the final col¬ 
umn. In the SED fitting for the MIR wavelengths, we place 
emphasis on the band fluxes (IRS-1,2,3,4,5) and flux densi¬ 
ties (IRAC-4 and MIPS-1) rather than the atomic line fluxes, 
because our interest in SED modeling is in calculating the 
gas and dust masses. Therefore, there are some discrepan¬ 
cies in the MIR atomic lines between the observed and calcu¬ 
lated data. The values are listed in the bottom line of Ta- 
ble[T4] The chi-square analysis implies that, within 1-cr, there 
was no difference between the predicted and the observed flux 
densities/band fluxes, but rather a slight (negligible) disagree¬ 
ment between the calculated and observed fluxes, owing to 
the C III] AA 1906/09 A flux. Without the C iii] AA 1906/09 A 
flux, x^ = 34.75 indicates that the modeled flux densities and 
band fluxes are in excellent agreement with the observations. 

The discrepancy between the observed calculated 
C III] AA 1906/09 A line fluxes appears to result from 
fluctuations in the structure of Tg. The Cm] lines are the 
most sensitive to the among those considered in the model; 
the excitation energy difference between the upper and the 
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Wavelength (yu,m) 

Figure 13. The SED of K648. The blue solid and broken lines are the predicted SED of the sum of the CSPN and the nebula and the incident SED of the 
CSPN, respectively. The black lines are the observed spectra obtained from the HST/FOS and Spitzer/IRS. The Subaru/HDS spectrum is not shown because 
its appeared to be partially contaminated by foreground stars. Instead, we plotted the HSTfW¥FC2 photometry, indicated by the green asterisks (CSPN) and 
triangles (CSPN+PN). The green triangles in the MIR and FIR (3.6/4.5/5.8/8.0/24/70/100/160 fim) are SpitzerflRCS/MlFS and HerschelfPACS photometry. 


lower levels (x) is 6.5 eV and the excitation temperature 
is 75 380 K {=x^k, where k is the Boltzmann constant). 
We used = 10270 K in the calculations of whereas 
the volume-averaged in the model was 11 090 K. 

With a constant n^, but a difference of only 820 K, the 
volume emissivity of this complex line at 11 090 K became 
^1.65 times larger than that with 10270 K. Accordingly, 
we obtained C Hi] AA 1906/09 A fluxes that were greater 
than those of the observations by a factor of ^1.65. The H/3 
emissivity at = 10270 K was 1.07 times greater than that 
at T): = 11 090 K. Therefore, the modeled C^'*' abundance 
was smaller than the observation by ^-0.10 dex. Taking the 
differences in the structure of Ug and between the model 
and the observed data into account, we estimate that the 
accuracy of the elemental abundances calculated using the 
was within ^0.2 dex. 

As the model provides decent predictions for the 
[Sill] AA 18.7/33.5Atm, [Siii]A9532A (/(P-I)) of 
[S III] A9532 A = 5.760, which are not listed in Tableand 
[S III] A 6312 A simultaneously, we may assume that the two 
MIR [S III] lines are not spurious, but rather genuine features 
of the spectra. 

Our P-I model with the CLOUDY code was not able to fit 
[Ne II] A 12.80 pm, whereas the prediction of the other atomic 
lines of similar IPs, i.e., ions such as [F ll] (see below), [S Hi], 
[Arm], and [Cl hi], is in good agreement with the observa¬ 
tions. Many P-I models using CLOUDY have been used to 
fit the [Ne ll] A 12.80 pm in PNe; however, to our knowledge , 
there has been little success (e.g.. lPottasch et ani201 lLl200^ . 
In our model, we monitored the chi-square values to obtain the 
best fitting parameters. With an almost constant radial density 


profile, the [Ne ll] A 12.80 pm could be modeled; however, the 
other line-fluxes and band fluxes/flux densities exceeded the 
observed values, i.e., chi-square increased. The recombina¬ 
tion rates for some heavy element ions (e.g.. S'*') are uncertain, 
so that photo-ionization models may give line fluxes that are 
in poor agreement with measured data. Therefore, the lack of 
agreement may be due to the uncertainties in the atomic data 
for Ne'*'. Improvements in these data, however, are beyond the 
scope of this paper. 

Our P-I model predictions provide good fits to two of the 
[F ll] line intensities. If both lines are not [F ll] lines but other 
elemental lines, the P-I model cannot fit these two lines si¬ 
multaneously. Therefore, we conclude that the detected [F ll] 
lines are likely to be real. The two observed [F ll] line fluxes 
and the calculated elemental abundance of F using the ICF(F) 
are in good agreement with the predictions of the model. 

The second and third columns of Table [16] list a compari¬ 
son of the nebular elemental abundances determine using the 
semi-empirical ICF method and the P-I model. As we men¬ 
tioned above, the accuracy of the elemental abundances de¬ 
termined using the P-I model was within ^0.2 dex. Careful 
treatment for the Ne abundance is necessary for the reasons 
discussed above. Therefore, we excluded the Ne abundance 
from the following discussion. The difference between the 
two data sets, (A), is listed in the fourth column of Table [Tbl 
The agreement between the He, C, O, F, S, Cl, and Ar abun¬ 
dances between the ICF method and the P-I model is generally 
good. 

However, poor agreement is found for N and Fe (|A|>0.2 
dex). The final two c olumns of Table [Tbl list the ICF values 
used in Section 13.1.51 and those predicted by the P-I model. 
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Table 15 

Comparison between the results of the P-I model and the observations. 


Ion 

A 

/(P-I) 

/(Obs) 

A///(Obs) 


(A//rm) 

[/(H/3) = 100] 

[/(H/3) = 100] 

(%) 

Cm] 

1906/09 

468.701 

334.984 

39.92 

ICii] 

2323 

25.134 

17.091 

47.06 

10II] 

3726 

17.646 

17.383 

1.51 

10II] 

3729 

9.385 

9.387 

0.02 

[Nelll] 

3869 

10.255 

9.939 

3.18 

[Nelll] 

3968 

3.091 

3.147 

1.78 

ClI 

4267 

0.492 

0.726 

32.23 

H 7 

4340 

46.974 

46.674 

0.64 

10 III] 

4363 

2.273 

2.782 

18.31 

He I 

4388 

0.640 

0.517 

23.84 

He I 

4471 

5.218 

4.914 

6.18 

[FII] 

4791 

0.105 

0.110 

4.94 

[FII] 

4870 

0.033 

0.030 

8.83 

[Fe III] 

4881 

0.054 

0.051 

4.92 

He I 

4922 

1.383 

1.290 

7.23 

[OlII] 

4931 

0.032 

0.038 

14.87 

[Olll] 

4959 

79.089 

74.974 

5.49 

[Olll] 

5007 

238.058 

227.263 

4.75 

[Cl III] 

5518 

0.022 

0.021 

3.14 

[Cl III] 

5538 

0.027 

0.028 

3.50 

[Nil] 

5755 

0.064 

0.052 

22.17 

He I 

5876 

15.676 

14.834 

5.67 

[Sill] 

6312 

0.147 

0.119 

23.55 

[Nil] 

6548 

0.979 

0.901 

8.70 

Ha 

6563 

282.047 

282.399 

0.12 

[Nil] 

6584 

2.890 

3.180 

9.12 

He I 

6678 

4.182 

4.114 

1.65 

[SII] 

6716 

0.069 

0.087 

20.98 

[SlI] 

6731 

0.110 

0.133 

17.23 

[Arm] 

7135 

0.389 

0.384 

1.21 

[On] 

7323 

1.573 

1.792 

12.20 

[On] 

7332 

1.256 

1.450 

13.39 

Hi 

7.47 

3.192 

3.150 

1.35 

[Arm] 

9.00 

0.296 

0.324 

8.64 

[SIV] 

10.51 

1.095 

1.064 

2.95 

[Nelll] 

15.55 

8.736 

11.545 

24.33 

[Ne II] 

12.80 

3.645 

14.980 

75.67 

[Sill] 

18.71 

2.343 

1.332 

75.93 

[Sill] 

33.47 

0.893 

0.612 

45.89 

IRS-1 

8.55 

13.365 

15.859 

15.73 

IRS-2 

9.825 

4.809 

4.560 

5.46 

IRS-3 

12.03 

4.973 

4.531 

9.75 

IRS-4 

14.00 

2.672 

2.322 

15.08 

IRS-5 

25.50 

9.044 

9.136 

1.01 

Band 

A 

/v(P-I) 

/v(Obs) 

/\FJFAOhs) 


(A//rm) 

(mJy) 

(mJy) 

(%) 

FI 60 BW 

1513 

26.647 

15.993 

66.61 

FI70W 

1820 

25.808 

20.886 

23.51 

F255W 

2599 

15.952 

11.907 

33.97 

F300W 

2989 

13.823 

10.543 

31.11 

F336W 

3360 

12.211 

13.393 

8.82 

F439W 

4312 

9.554 

9.793 

2.43 

F547M 

5484 

5.810 

6.025 

3.56 

F814W 

7996 

3.701 

3.921 

5.62 

IRAC-I 

3.51 

1.273 

5.096 

75.01 

IRAC-2 

4.50 

1.499 

4.158 

63.95 

IRAC-3 

5.63 

2.035 

5.046 

59.67 

IRAC-4 

7.59 

4.734 

8.510 

44.37 

MIPS-1 

23.21 

11.300 

10.684 

5.77 

PACS-B 

68.93 

3.207 

2.950 

8.71 

PACS-R 

153.9 

1.804 

2.680 

32.69 





88.12 


Note. — The data in the IRS-1, 2, 3, 4, and 5 bands ai‘e the integrated 
fluxes between the following wavelengths: 8.26-8.84 /im, 9.7-9.95 ^m, 11.9- 
12.16/im, 13.9-14.1/im and 24.5-26.5 /im, respectively. Data are shown 
with two or three decimal places to avoid rounding eiTors. 


The ICF values from the P-I model were generally in agree¬ 
ment with those of the semi-empirical methods, except for N. 
This is because most fractional ionizations occurred in other 
ionic stages: the P-I prediction suggests 5% for and 95 % 
for Poor agreement for N between the ICF and P-I mod¬ 
els is oft en found for PNe and in the O-rich halo PN DdD ml 
(see e.g.. lOtsuka et al.|[2009l : lDelgado-Inglada et al]l2014f) . In 
many cases, including DdDml, N'*' abundances determined 
from optical spectra alone have been used to determine the 
elemental N abundance, because it is difficult to detect the 
forbidden lines, which appear in t he UV or FIR spectra. 
Based on grid models using CLOUDY. iDelgado-Inglada et al.l 
(120141) proposed that the N/O ratio for PNe showing no-He II 
lines can be estimated by the following equations, 

N N+ 

-=ICF(N)Gn4-^, (5) 

ICF(N)gii 4 = 10° o+°+o2+ . (6) 

In the case of K648, the ICF(N)gii 4 and the N/H abun¬ 
dance using the observed abundances (Table |7]i 

and the elemental O abundance (Table [TOb were estimated to 
be 3.40±2.42 and 7.77(-6)±5.82(-6). The logio(N/H)H-12 
of 6.89±0.33 is very close to the model predicted value 
(6.96±0.20), although we should note that the ionic and el¬ 
emental abundances would depend on the density structure, 
incident ionization source, ionization boundary condition, gas 
metallicity, dust grains/molecules, and so on. We should keep 
in our mind that we were unable to detect N lll|^A1750 A in 
K648. The predicted line-intensity around 1750 A FOS spec¬ 
trum with S/N = 1 was ~5.3 (i.e., a detection limit), which is 
approximately three times larger than the prediction of the P-I 
model (i.e., 1.7). The N III] A 1750 A line was too faint to de¬ 
tect using FOS. As the P-I model provided a value of ICF(N) 
that was too large, in this paper we prefer to use the semi- 
empirically determined N abundance, i.e., the ICF abundance 
(not ICF(N)gii 4 ). The P-I model indicates that Fe ions are also 
concentrated in other ionic stages, rather than the observed 
ionic stage, i.e., 20% Fe^"^ and 80% Fe^"^. The availability of 
the and Fe^^ lines would be expected to improve the ac¬ 
curacy of the abundance calculation. FIR observations using 
^P/CA/SAFARI would be helpful to detect [N ll] A 121.3 /rm 
and [N III] A57.3 /rm lines and verify ICF(N) in PNe. 

Note that we used the ICF abundances rather than the P-I 
results. Analysis with the P-I results, however, is not expected 
to alter the conclusions. The P-I results should be carefully 
examined in a more sophisticated future study. 

The nebula is full y ionized, so that the ionized gas mass is 
consistent with mg. iBianchi et alJ (119951) determined an ion¬ 
ized g as mass of 0.05-0.09 Mq, and iKingsburgh & Barlow! 
(119921) determined it to be 0.042 Mq. Both authors assumed 
a constant density profile. Although mg (and mj) depends on 
the distance used, mg in this work is consistent with these es¬ 
timates. 

Here, we estimate and the dust-to-gas mass ratio md/nig. 
If we assume that the AC dust grains have a radius of 0.0005- 
0.25 fim and a size distribution that follows we obtain 
mj = 9.74(-7) Mq and moling = 2.07(-5). However, this 
model does not fit the flux density at the above-mentioned 
wavelengths. To our knowledge, the value of nid = 4.95(- 
7) Mq for K648 is the smallest mass among known PNe, 
and is approximately one order of magnitude smaller than 
that of BoBnl, where the nij = 5.78(-6) Mq. For BoBnl, 
lOtsuka et al.l (120101) used grains with a radius of 0.001- 
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Table 16 

Comparison of the observed elemental abundances and those predicted using the P-I model. 


X 

Obs“ 

log(X/H)+12 

P-P 

log(X/H)+12 

log{Xobs/Xp-i) 

ICF(XobA^ 

ICF(Xp_i)= 

He 

11.02±0.03 

10.99±0.20 

-b0.03±0.20 

1.04±0.01 

1.00 

C 

8.97±0.17 

8.71±0.20 

-b0.26±0.26 

1.33±0.24 

1.00 

N 

6.36±0.10 

6.96±0.20 

A).60±0.22 

4.02±0.81 

20.56 

O 

7.73±0.03 

7.82±0.20 

-0.09±0.20 

1.00 

1.00 

F 

5.42±0.11 

5.41±0.20 

+0.01 ±0.23 

4.02±0.81 

5.22 

Ne 

7.44±0.03 

7.02±0.20 

+0.42±0.20 

1.00 

1.00 

S 

5.40±0.07 

5.48±0.20 

A).08±0.21 

1.00 

1.00 

Cl 

3.58±0.15 

3.44±0.20 

+0.14±0.25 

1.19±0.29 

1.05 

Ar 

4.60±0.13 

4.44±0.20 

+0.16T0.24 

1.19T0.29 

1.04 

Fe 

5.02±0.12 

5.48±0.20 

-0.46±0.23 

4.02±0.81 

5.15 


^ From Table|T^ We used the RL He and the CEL C/N/O/F/Ne/S/Ar/Cl/Fe abundances in the P-I model. 
^ Determined from the P-I model. 

^ Elemental abundance difference between the observed and the model predicted abundances. 

^ from TablelTol 
^ calculated from the P-I model. 


Table 17 

Expansion velocities of K648. 


Ion 

Type of 
lines 

IP. 

(eV) 

Num. of 
sample lines 

(Vexp) 

[Oi] 

CEL 

0.00 

2 

10.12±0.27 

[SII] 

CEL 

10.36 

2 

14.70±0.28 

Hi 

RL 

13.59 

26 

15.07±0.58 

[On] 

CEL 

13.62 

2 

15.99±0.18 

[Nil] 

CEL 

14.53 

3 

13.81±0.43 

[Fe III] 

CEL 

16.18 

2 

12.70±1.05 

[Eli] 

CEL 

17.42 

2 

13.85±0.90 

[Sill] 

CEL 

23.33 

1 

15.96T0.57 

[Cl III] 

CEL 

23.81 

2 

16.12±1.55 

CII 

RL 

24.38 

6 

19.37±0.99 

He I 

RL 

24.59 

17 

15.71±0.56 

[Arm] 

CEL 

27.63 

1 

14.25±0.16 

N II 

RL 

29.60 

2 

15.72±1.49 

[Olll] 

CEL 

35.12 

4 

15.75T0.52 

Oil 

RL 

35.12 

1 

12.65T0.99 

[Nelll] 

CEL 

40.96 

2 

13.75±0.09 

Cm 

RL 

47.89 

3 

15.54±0.63 


Note. — The third and fourth columns list the IP and the number of sample 
lines used in the calculation of the average Eexp of each ion, respectively. 


0.25 Jim with an a~^-^ size distribution in the SED model, 
whereas for K648 we used much smaller grains to match the 
observed SED at wavelengths in the range 10-15/rm. We 
found the ratio nijlnig = 1.03(-5), which was much lower 
than th at for BoBnl (5.84x10“^). ForH4-l. lTaiitsu & Otsukal 
(l2014h reported that mg = 0.3 Mq, = 7.34(^) Mq, and 
mjlnig = 2.48(-3); however, H4-1 contains abundant cold dust 
and hydrogen-rich molecules. Among these C-rich halo PNe, 
where the metallicity is similar to that of K648, we could not 
find dependence of the metallicity on the dust mass. 

3.5. Expansion velocities and the time since the AGB phase 

We employed a multiple Gaussian fitting method for 
the flux measurements, except for the strong lines 
[Oh] AA 3726/29 A, [Oiii] A A 4959/5007 A, [Oi]A6300A, 
Ha, and [N ll] A 6583 A, because these lines have a weak 
broad tail component or a small offset velocity component. 
Here, we focus on the nebula expansion velocity Vexp deter¬ 
mined from the main Gaussian component of each line. 

We measured Vexp using the following relation; 

Vexp = 1/2 (VFWHM“V,heim“ (7) 

where Vfwhm is the FWHM of the velocity, Vtherm is the ther¬ 


mal bro adening velocity, and Vns t r is the instrumental veloc¬ 
ity (e.g.. lOtsuka et al.ll2010ll200^lBianchi et al.ll2001l) . Here, 
Vtherm is represented by 2\A{T^'K\(r^/Ar)^!^, where is the 

relative atomic mass of the target ion. For CELs, we used 
the Tg listed in Table |6] For RLs, we used for H I, 

CII, III, [Nil], and Oil and the Tg/Hel) = 6710 K for the 
He I lines. We measured Vinstr for all the identified lines 
listed in Table in the Appendix using the nearby Th-Ar 
lines, i.e., 4.3 km s”'for [Arm] A7135 A, Hei A7281 A, and 
[OII] AA 7320/7330 A (the resolving power of these lines was 
^69000), and 8.8-9.0 km s“' for the others . We did not 
include the turbulent velocity, because these vel ocities have 
been m eas ured in ^^100 Galac tic PNe by e.g., lAcker et"^ 
and iGesicki et akl (l2003h . who found no turbulent ve¬ 
locities in PNe with non-WC type central stars, such as K648. 

The resulting Vexp are summarized in Table [IT] We mea¬ 
sured the VexD of over 100 lines selected from the lines listed 
in Appendix TablelAI For each ion, we excluded the measure¬ 
ments far from the average using l-cr clipping. We then cal¬ 
culated the average expansion velocity, (Vexp) of, the 18 ions 
listed in the final column of Table [TtI 

Our measu r ement s showed good agreement with those of 
iBianchi et al.l (1200Ih . who measured the (Vexp)(Hl) of 16.7 
km s“* using Ha and (Vexp)([N ll]) of 11.9 km s“*, using both 
[N II] A 6548 A and [N ll] A 6583 A lines, and with a constant 
electron temperature of Tg = 10000 K for the thermal broad¬ 
ening velocities. We used 26 Hi lines for the (Vexp)(Hl) 
and 3 lines for the (V'exp)([N ll]) calculations with 7’g(BJ) and 
7’g([NH]), respectively. Th e slight differences b etween our 
data and those reported by IBianchi et alJ (1200 Ih can be at¬ 
tributed to the value of Tg used and the number of sample 
lines. 

The correlation between the (Vexp) and IPs is given by 

(Vexp) =(5.23 ± 4.03)-10”^ IP3-13.54 ± 1.01. (8) 

The correlation factor was 0.32. Assuming that K648 has 
a standard ionized structure, i.e., high-intensity IP lines are 
emitted from regions close to the central star and low-intensity 
IP lines are emitted from regions far from the central star, the 
expansion velocity of the nebula may be slowing with an al¬ 
most constant value of r. In general, the PN shell is known to 
follow a Hubble type expansion, i.e., acceleration of the ex¬ 
panding gas shell. Perhaps it did not gain its impulsion from 
the CSPN yet. 

As we found for BoBnl (lOtsuka et alJl201(il) . Vexp for [O ll] 
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was at least 1.5 km s“^ smaller than that for [O III]. 

The apparent outer radius of K648 is 2.1", which corre¬ 
sponds to 0.125 pc at 10.9 kpc. The wind velocity of K648 
during the AGB mass-loss phase is unknown. We used an 
expansion velocity of (Vexp)(Hl) = 15.07 km s“', and esti¬ 
mated the dynamical age of the K648 nebula to be 8110±490 
years since the AGB phase. iRauch et"^ (120021) estimated the 
post-AGB age of 6800 ^|qqq yrs by plotting their derived lu¬ 
minosity and surface gravity on theoretical evolutional tracks. 
However, the evolutionary age after th e AGB phase with more 
precision is unknown at this moment. [McCarthy et alJ (119901) 
investigated the disagreement between evolutionary and dy¬ 
namical time scales for the evolution of the CSPNe using 
the results of high resolution spectra of about 23 CSPNe. 
According to them, the AGB-CSPN evolutionary transition 
times could have been increased by small additional amounts 
of residual envelope material remaining after the superwind 
mass-loss phase. 

4. DISCUSSION 

4.1. Comparison with the AGB nucleosynthesis model 

In Section 13.2.31 we determined the core-mass of the cen¬ 
tral star as 0.61-0.68 Mq, depending on the choice of the dis¬ 
tance to M15. The initial-final mass relation has been stud¬ 
ied using solar metallicity fo r young (^^1-2 Gyr) open clus¬ 
ters (e.g.. lKalrrai et alJl2008l) : however, it has been not stud¬ 
ied using metal-poor old clusters. Semi-empirical initial-final 
mass relations are only available for the chemical compo- 
sitio n of the solar neighborhood and for Magellanic Clouds 
(see. lPrada Moroni & Stranieroll200'^ . The mass-loss and the 
dredge-up efficiency (depending on the core-mass, metallic¬ 
ity, and total mass of the star) during the AGB phase de¬ 
termine the fate of stars. From these reasons, we utilized 
the theoretical initial-final ma s ses for Z = 10“^ stars reported 
by iPrada Moroni & Stranierol (120071) to estimate the initial 
mass of K648. From pol ynomial fitting to the init i al-fina l 
masses listed in Table 1 of iPrada Moroni & Stranierol (l2007h . 
we found that core-masses of 0.61, 0.63, 0.66 and 0.68 Mq 
correspond to the initial masses of 1.15, 1.60, 1.76, and 1.87 
Mq, respectively. As the upper limit of the mass of stars in 
M15 is ^1.6 Mq, the current core-mass and initial mass of 
K648 would be ~0.61-0.63M^and ~1.15-1.6 Mq, based on 
iPrada Moroni & Stranierol (1200^ . 

For comparison, we discuss the results of iLugaro et alJ 
(120121) for 0.9, 1.25, and 1.5 Mq stars with an initial [Fe/H] 
= -2.19. ILugaro et al.l (120121) used scaled solar abundances 
as the initial composition for all elements from Li to Pb 
[X/Fe]~0 and [He/Fe] = -(-2.18. The initial conditions and 
the mass l oss formulae us ed in ILugaro et alJ (12012^ were dis¬ 
cussed bv lKarak^ (1201 Oh . Table ITS] lists the predicted abun¬ 
dances after the final TP. Here, we used the nebular abun¬ 
dances, except for P, where we used the stellar abundance. 
The C and O abundances used for K648 were the values from 
the CELs. The 0.9, 1.25, and 1.5 Mq stars would, theoreti¬ 
cally, experience 38, 15 and 18 TPs, respectively. The final 
three lines of Table fTSl lists the ejected mass during the final 
TP, the final core-mass, and the envelope mass. 

Our estimated core-mass agrees with the pre dictions for 
1.25 and 1.5 Mq reported bv lLuigaro et al.l (120121) . As the 0.9 
Mq models experienced many TPs, the final core-mass was 
larger than that pre dicted by the 1.25 and 1.5 Mq models. 
ILugaro et alJ (120121) included a partial mixing zone (PMZ), 
which is formed in a mixing zone from the H-rich enve¬ 


lope down to the layer at the top of the He-rich intershell 
(IShingles & KarakM 120131) . The PMZ produces a (as 
well as a ''^N) pocket during the interpulse period. The 
releases additional free neutrons (n) via ^^C(a,u)'®0, result¬ 
ing in further «-process elements, such as '^F and ^'P. The 
available mass of mainly affects the final compositions of 
K648, in particular, C, N, O, Ne, and F, which are syn thesized 
in the He-rich intershell. IShingles & Karakai (120131) showed 
that the Ne abundance increases as the mass of the PMZ in¬ 
creases. They argued that the Ne enhancement is due to ^^Ne 
production via double a-particle capture by '^N. 

For the reason given above, we checked the abundances of 
N and Ne. First, the 0.9 Mq and 1.25 Mq model with no PMZ 
were excluded. The former could explain the Ne abundance, 
but not the N abundance. In addition, the final core-mass ap¬ 
peared larger than the predictions of these models. The latter 
model could not explain the Ne abundance either. The re¬ 
maining two, i.e., the 1.25 Mq+2(-3) Mq PMZ and 1.5 Mq 
models, provided reasonable agreement with not only the N 
and Ne abundances, but also the He, C, O, and F abundances. 
The P abundance of the CSPN was comparable to the pre¬ 
dicted value of 3.49 using this model. 

Our estimate of the core-mass of the CSPN is in good agree¬ 
ment with both the 1.25 Mq and 1.5 Mq models; however, 
the resulting gas masses could not be explained using either 
model. Therefore, we expect that models for stars with initial 
masses in the range of 1.25-1.5 Mq can explain the ejected 
mass, as well as elemental abundances and the core-mass of 
the CSPN. 

4.2. Was the progenitor a blue straggler? 

We found that even the lower core-mass of 0.61 Mq, which 
corresponds to an initial mass of 1.15 Mq, exceeds the mass 
of turn-off stars in M15; however, the uncertainty of ^0.03 
Mq should be noted. The 1.25 Mq+2{-3) Mq PMZ model 
can also explain observed nebular abundances. Hence, it is 
possible that the progenitor of K648 is a binary system. In¬ 
deed, K64 8 has long been sus pected to have undergone binary 
evolution (I Jacoby et al.lll99'^ . 

During the evolution of the progenitor of K648, if it effi¬ 
ciently gained mass and nucleosynthesized products via mass- 
transfer from an evolved massive primary, it could evolve into 
a C-rich PN. Although it initially appears difficult to accept 
that binarity may be responsible for many cases of anoma¬ 
lous composition, there is now evidence of radial velocity 
variations or bright equatorial disk structures, signaling a bi¬ 
nary orbit, and binary interactions are regarded as the ex¬ 
planation of a wide range of C-rich stellar classes, barium 
stars, CH stars, an d CEMP stars (Carbon-Enhanced Metal- 


Poor stars, see, e.g.. 

Beers & Christliebll2005tlMasseron el alJ 

l2010l;lBisterzo et al. 

[ToiHv 


In th e case of the C-rich halo PNe H4-L lOtsuka & Taiitsul 
(120131) proposed that the H4-1 chemistry may be the evolu¬ 
tionary result of a ~0.8-0.9 Mq star that had been affected by 
mass-transfer from a more massive AGB comp anion in a bi¬ 
nary system; however, lOtsuka & Taiitsul (120131) did not eval¬ 
uate the core-mass of the central star due to the lack of obser¬ 
vation data of the central star. The traditional evolution theory 
that the progenitors of PNe are the remnants of single stars at 
the end of the AGB phase does not provide a natural expla¬ 
nation for the non-spherical morphologies observed for the 
great majority of PNe. Although the binary interaction model 
explains some of the anomalies associated with the observed 
PN population, the number of PN central stars with known 
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Table 18 

Comparison of the observed nebular abundances with the predictions of the [Fe/H] = -2.19 AGB model. 





Models 




Obs 


Initial mass (Mq) 

0.9 

0.9 

1.25 

1.25 

1.5 


Elements 

PMZ mass (M©) 

0 

2(-3) 

0 

2(-3) 

0 


He 


11.00 

11.00 

11.00 

11.00 

11.01 

11.02±0.03 

C 


9.07 

9.04 

8.94 

8.90 

9.26 

8.97±0.17 

N 


7.55 

7.53 

6.67 

6.68 

6.76 

6.36±0.10 

0 


7.48 

7.63 

7.31 

7.47 

7.56 

7.73±0.03 

F 


4.74 

5.03 

4.35 

4.78 

5.08 

5.42±0.1I 

Ne 


7.33 

7.87 

6.95 

7.65 

7.68 

7.44±0.03 

pa 


3.48 

3.64 

3.36 

3.54 

3.49 

3.64±0.10 

ejected mass during 


2.0(-3) 

2.0(-3) 

8.0(-3) 

8.0(-3) 

5.99(-l) 

4.8(-2)‘’ 

last TP (M©) 








core-mass (Mq) 


0.77 

0.77 

0.66 

0.66 

0.66 

0.61-0.63 

envelope mass (Mq) 


0.04 

0.04 

0.02 

0.03 

0.18 



Note. — The initial He, C, N, O, F, Ne and P abundances in all models are 10.92, 6.29, 5.69, 6.55, 2.24, 5.79 and 3.23, respectively. 
“ The P abundance of the CSPN measured in the FUSE spectrum. 

The mass estimated using the CLOUDY SED model. 


binary companions is very small and carrying out programs 
to det ect such objects are extremely difficult (see Llacobv et alJ 
120131 and references therein). 

H4-1 and K648 both appeared to exhibit evidence of binary 
evolution structures, s uch as bright equatorial dis k structures 
and a bipolar nebula (iTaiitsu & Otsukal 120041) . lAlves et aP 
(120001) did not detect any time variation in the magnitude of 
the central star using //57’AVFPC2. They argued that the fail¬ 
ure to detect a current binary companion lends support to the 
picture of a complete merger, as opposed to more modest 
mass transfer, because in the latter case there would still be 
a remnant companion (possibly a helium-rich white dwarf). 
We neither detected any variation in the radial velocity in our 
HDS spectra. Although the failure to detect a companion is 
not conclusive proof, it is worthwhile to re-examine whether 
the central star of K648 is a binary or a merger. 

K648 might be not a “typical” PN. For example, the ion¬ 
ized gas mass is unusually small, which may indicate that it 
formed via a non-standard mechanism. Based on a compari¬ 
son with AGB yields of a single star, as discussed above, we 
propose that the merging of two stellar bodies occurred, or a 
large mass fraction was transferred from a companion. In a 
close binary system, the gravity of one component can induce 
a significant tidal force in the other. The dissipation of this 
tidal force may synchronize the rotation and circularize the 
orbit, leading to coalescence (or consumption the outer enve¬ 
lope of its companion) in extreme cases. Although there have 
been many theoretical analyses and simulations of binary co¬ 
alescence of neutron stars or black holes, there have been no 
reports of clos ely related work in bin ary systems in PNe such 
as K648 (e.g.- IZhang & .Teffervll2013l and references therein). 

The notion that K648 and other globular clusters may 
arise from c oalesc ence of binary syst e ms w as proposed by 
lAlves et ^ (120001) and iJacobv et al.l (Il997h . lAlves et al] 
(l2000h at^ed that the progenitor of K648 experienced mass 
augmentation in a close binary merger, and evolved as a 
higher mass star to become a PN. Such a high-mass star would 
be a blue straggler (BS). A number of pos sible BS candidates 
(20-69 objects) have been found in Ml5 (iDieball et al.ll200^ 
iDiaz-Sanchez et ^l2012h . 

There are several ways in which stars may evolve into BSs, 
i.e., MS-MS collisions, WD-MS collisions, and close binary 
transf er or mergers (e.g., llJmbreit et alJ I200M iFerraro et alJ 
l2009h . The progenitor of K648 may have been formed via 


a close orbital activity of a binary with a large mass inflow 
from its companion during the MS stage, one possibility is 
a close binary system that consists of two stars of Q.9Mq 
with slightly different masses, or with sign ificantly different 
masses as proposed bv lOtsuka et alJ (1201 Oh . and that one star 
consumed a large mass fraction of its companion so that the 
mass of this star would approach ^1.6 Mq, relegating its 
companion to the position of an accessory. 

To date, 20-69 BSs (including candidates) have been 
identified in M15 (IDieball et al.l 120071 : iDfaz-Sanchez et aP 
2012h. The typical PNe lifetime is -^ 25 000 years (e.g., 
Moe & De Marcoll200^lFel dmeieil2 (ffl3l) . When we use a BS 
lifetime of ~ 1.2 Gvrs dsills et al.ll 19971) . the expected number 
of PNe in M15 is 0.0004-0.0014 PNe per a BS (=25 000 years 
/ 1.2 Gyrs x 20-69). If we assume a birth rate of 2.5-5.0(-8) 
BSs yr“^ via this process (lUmbreit et alJl20()8l) . and given that 
the age of M15 is 13.5 Gyrs, the estimated number of PNe 
formed in M15 is 0.135-0.675 PNe. Therefore, K648 would 
be a rare PN evolved from a BS. The central star of K648 
may be a BS of higher mass determined by us (i.e., 1.5 Mq) 
or close to the limit in M15 (^1.6 Mq). Other BSs found 
in M15 may well evolve into PNe, similar to K648. If new 
evidence of much lower stellar mass is reported, the milder 
binary interaction scenario must be explored accordingly, i.e., 
mass-transfer from a more massive AGB companion in a bi¬ 
nary system. Flowever, this scenario may be not appropriate 
to explain the relatively large mass of K648. 

4.3. Comparison ofK648 with BoBnl and H4-1 

lOtsuka et alJ (1201 Oh reported a similar comparative study 
of their data for BoBnl with earlier analyses of K648. Based 
on inco mplete observations, e.g., no detection of [Fll] lines 
in K648. IOtsuka et aP (1201 Oh argued that BoBnl might have 
undergone binary evolution with a 0.75 Mq+\.5 Mq system, 
whereas K648 might be an object that went through either a 
binary evolution with 0.75 M0-(-1.5 Mq or a single 1.8 Mq 
stellar evolution, ignoring the upper limit of mass for stars in 
M15. Here, we refined the earlier guessing based on the abun¬ 
dances of 10 elements. In contrast to BoBnl, no evidence of 
Ba and Xe was observed in K648; however, there are simi¬ 
larities between BoBnl and K648 which are not shared with 
H4-1. Here, we discuss similarities between K648, BoBnl 
andH4-l. 

Figure[T4]shows the elemental abundance patterns of K648, 
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Elements (X) 

Figure 14. The relative elemental abundances of K648, H4-1 and BoBnl 
com pared with those of the Su n. T he solar Kr, Xe, and Ba abundances are 
from ILoddersI <200-3l> . See Table ll9l for the elemental abundances of each PN. 


Table 19 

Comparison of the abundances for K648, H4-1, and BoBnl. 


Elements 

K648 

H4-1 

BoBnl 

Average 

He 

1L02±0.03 

1L03±0.15 

1L07±0.01 

1L04±0.06 

C 

8.97±0.17 

9.02±0.18 

9.02±0.08 

9.00±0.14 

N 

6.36±0.10 

7.59±0.04 

8.03±0.10 

7.69±0.09 

O 

7.73±0.03 

8.18±0.02 

7.74±0.03 

7.94±0.02 

F 

5.42±0.03 


5.85±0.09 

5.68±0.09 

Ne 

7.44±0.03 

6.43±0.10 

7.96±0.02 

7.60±0.03 

S 

5.40±0.07 

5.13±0.03 

5.32±0.16 

5.30±0.09 

Cl 

3.58±0.15 

3.88±0.13 

3.39±0.07 

3.6640.12 

Ar 

4.60±0.13 

4.56±0.12 

4.33±0.04 

4.5140.10 

Fe 

5.02±0.12 


5.08±0.13 

5.0540.13 

Kr 



2.88 

2.88 

Xe 


>2.70 

<2.97 

2.86 

Ba 


<2.51 

1.97 

2.32 

C7D 

I7.47±7.07 

6.93±2.96 

19.06±3.75 

14.4944.59 


Note. — The abundances of all elements except He were determined from 
the CEL lines; that of He was determined from RL lines. The elemental abun¬ 
dances are in the form of log[o(X/H)-l-12, where H is 12. The C/0 number 
density ratio is linear va lue. The values of H4-1 and BoBnl were taken from 
lOtsuka & Taiitsul <201311 and lOtsuka et all <201(11 . 

as well as those of H4-1 and BoBnl. The abundances of all 
elements except He were determined based on CEL lines; the 
abundance of He was determi ned from the RL l ines. Data 
for H4-1 and BoBn l are from lOtsuka & Taiitsul (1201 3h and 
lOtsuka et al.l (1201 Ol) . respectively. Discrepancies in the C and 
O abundances in H4-1 and BoBnl are discussed in these re¬ 
ports. A comparison of the RL C and O abundances among 
K648, H4-1 and BoBnl is beyond the scope of this paper. The 
elemental abundances for each PN and the average abundance 
of each element are listed in Table [19] 

The He and C abundances were the same for the three sys¬ 
tems to within error. The a-elements Ar and S, and Cl were 
not synthesized in significant quantities in the PN progeni- 
to rs. Lor examp le, with the 1.25 Mq -(-2.0(-3) PMZ model 
of iLugaro et al.l (120121) . an increase of only ^0.02-0.03 dex 
was found compared with the initial abundances. Therefore, 
we can regard these three elements as mostly SN products. 
Some fractions of O and Ne are synthesized in the He-rich 
intershell during the TP-AGB p hase. Indeed, the above 1.25 
Mq model iLugaro et al.l (120121) predicted that such stars can 
increase +0.94 and 4-1.88 dex from the initial O and Ne abun¬ 
dances, respectively. 


We determined gas-phase abundance of [Pe/H]=- 
2.45±0.12 using the ICP(Pe) in K648, which i s comparable 
to the typical [Le/H] abundance in a M15 star (ISobeck et alJ 
[2(iril [Le/H]—2 .3). There are a number of other possible 
forms of Le in the solid phase in laboratory experiments. 
According to the Jena Database of Optical Constant^ for 
example, LeO, LeS, m agnesium-iron oxides, magnesium-iron 
silicates, and olivine. iDelgado-Inglada & RodrfguezI (1201 4t) 
employed a more sophisticated ionization correction factor 
scheme than that used here. They reported that the highest 
depletion factors are found in C-rich objects, exhibiting SiC 
around 11 pm or a broad 30 pm feature in the infrared spectra. 
Note that the carriers of these features are under debates. 
The central positions of the S i C feat ures in the sample of 
iDelgado-Inglada & RodrfguezI (120141) were not reported, 
however, and K648, H4-1, an d BoBnl were not included in 
their s ample. According to IDelgado-Inglada & RodrfguezI 
([20I3), the Le abundance ratio depletion detected in most 
PNe might be due to situation that less than 10 % of the Le 
is in the gas phase with more than 90 % in the solid phase. 
Instead of using the semi-empirical ICL scheme, we have 
an alternative method of estimating the [Le/H] abundance, 
i.e., using the CLOUDY P-I model. The [Le/H] abundance 
predicted using our P-I model was [Pe/H]=-1.99±0.2, which 
is close to that determined using the ICL method. The 
agreement was significantly b etter than that reported by 
IDelgado-Inglada & RodrfguezI ((MQ)- The resulting ICL 
[Le/H] value for K648 is likely to represent the abundance 
of Le, because this value was consistent with that for typical 
Ml5 metallicity, and we did not find features corresponding 
to amorphous silicate, crystalline silicates, or SiC (see section 

IlOi . 

There have been no reports of the detection of these features 
in H4-1, BoBnl, or K648. The presence of MgS and LeS is 
known to result in broad features around 30 pm (although the 
carrier of this feature remains the subject of some debate). 
We did not detect any other refractory element lines, such 
as Si and Mg, to estimate their ionized abunda nces in these 
PNe. The 1.25 MqH-PMZ 2(-3) Mq model of lLugaroet^ 
(1201 2^ predicted that the [Si/H] and [Mg/H] abundances are - 
2.15 and -1.67, respectively. If some of the Si and Mg-atoms 
might exist as dust grains, the gas-phase abundances of these 
two elements would become smaller, so it would be difficult 
to detect ionized emission-lines of these elements. The large 
fraction of Le and S cannot be due to dust grains such as MgS 
and LeS. Therefore, the S and Le abundances represent S and 
Le in these halo PNe, and these elements are expected to exist 
mostly in the gas phase. However, we cannot completely ex¬ 
clude the possibility that some fraction of Le resides in other 
solid forms. 

The abundances of S, Ar, Cl, and Le for K648 are approxi¬ 
mately the same as those for H4-1 and BoBnl. At first, it nay 
be expected that all three progenitors were born in the same 
chemical environment during the same epoch. However, there 
appear to be subtle differences in the birth environments as 
well as the evolutionary histories. 

The enhanced abundances of O, Ne, and the u-capture el¬ 
ements provide clues regarding the chemical environment 
where the progenitors originated and the nucleosynthesis in 
the inner core of the progenitors. Intrinsically larger [O/Le] 
appears in metal-poor stars, which is known to be the result 
of the time delay effects. However, note that the observed O 

* littp://www.astro.uni-jena.de/Laboratory/Database/jpdoc/index.html 
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abundance in our sample is the sum of SN and AGB nucle- 
osynthesized values. The O abundances in K648 and BoBnl 
are approxi mately equal. For H4-T however, after a detail 
discussion, lOtsuka & Taiitsul (1201 3h argued that 0.2-0.3 dex 
of the observed a-elements are SN products. The C/O ratio 
of K648 (the C-richness indicator) is very similar to that of 
BoBnl. The similarities of these two elements in both PNe 
may be explained as binary evolution and chemical enrich¬ 
ment during the AGB phase. Due to the O-richness, the C/O 
ratio of H4-1 is lower than that for the other two PNe. 

The Ne abundances vary significantly amongst these halo 
PNe. The very small Ne abundance for H4-1 indicates that 
the progenitor of this PN has no PMZ. A PMZ may have 
been formed in the He-rich intershells of K648 and BoBnl. 
For BoBnl, the Ne enhancement would be due to '"^N in the 
large PMZ and in H-burning ashes, and the ^^Ne enhancement 
via double a-capturing by '"^N. The enhancements of Ne and 
N in BoBnl are similar ([N/H] = h- 0.20±0.11 and [Ne/H] = 
-i-0.09±0.10), although AGB models do not yet reproduce the 
N and F overabundances in BoBnl. Nonetheless, the models 
successfully explain both the N and F abundances in K648. 
K648 therefore is expected to have been born in an similar 
chemical environment as BoBnl, but the progenitors experi- 
e nced different nucleosynt hesis. 

lOtsuka & Taiitsul (1201 3h argued that the abundance of Xe 
in H4-1 appears heavily polluted due to the r-process in Type 
II SNe, whereas the Xe abundance in BoBnl is close to the 
theoretically predicte d amount via 5-proce sses in AGB nucle¬ 
osynthesis; therefore. lOtsuka et al.l (1201 Oh concluded that the 
Xe in BoBnl is a product of s-processes. Therefore, the chem¬ 
ical environments where H4-1 and BoBnl were formed were 
very different. 

5. SUMMARY 

We have described observations of the PN K648 in M15 and 
investigated chemical abundances in the nebula, the CSPN, 
and dust-based regions, using multiwavelength data. We de¬ 
termined 10 elemental abundances for the nebula, including 
those for F, Cl, and Fe, which are reported here for K648 for 
the first time. The F enhancement in K648 is comparable to 
that for the C-rich halo PN BoBnl. We determined the C and 
O abundances from both CELs and RLs. The RL C abundance 
was consistent with the CEL value, whereas the RL O abun¬ 
dance was approximately three times larger than that of the O 
CELs. We attempted to obtain Ne abundance more accurately 
by adding the Ne"^ abundance determined using the Spitzer 
data. We determined the abundances of He, C, N, O, Ne, P 
and Ee, as well as the physical parameters of the CSPN, by 
employing a spectral synthesis fitting method. We found that 
the C/O and Ne/O ratios of the CSPN are roughly consistent 
with those of the nebula determined from the CELs and RLs 
within the error. The similar C/O ratios might indicate that 
the nebular abundances are reflective of the most recent stel¬ 
lar wind ejection from the central stellar surface. Spitzer/IRS 
shows the Class B 6-9 pm and 11.3 pm PAHs, as well as the 
broad 11 pm feature in K648. 

We constructed a CLOUDY radiative transfer P-I model to 
investigate physical conditions of the gas and dust in a self- 
consistent manner, and estimated the respective masses. The 
observed chemical abundances and core-mass of K648 are in 
agreement with AGB nucleosynthesis models for initial 1.25 
MqH-PMZ = 2x10“^ Mq stars, as well as initial 1.5 Mq stars 
without PMZ. Our simulation result confirms a possibility that 
K648 had evolved from a star with a mass in the range of 1.25- 


1.5 Mq. Perhaps the progenitor of K648 experienced coales¬ 
cence (or a large mass-transfer from its companion) during 
the early stages of evolution, and became a 1.25-1.5 Mq 
blue straggler (BS). If K648 is a PN that evolved from such a 
BS in M15, it would be a very rare or the first such case iden¬ 
tified among BS stars in Ml5, given that the expected number 
of PNe that evolved from BSs to date is only 0.135-0.675. 

We performed the analysis of all observational data avail¬ 
able, across a wide range of instruments and telescopes from 
the UV to the infrared for K648 with the help of P-I model 
construction. Based on our analysis, we proposed that K648 
could be evolved from a BS. The BS evolution scenario into a 
C-rich PN is still at the speculative stage. A detailed hydrody¬ 
namic simulation may help to visualize the population-based 
chemical evolution, or assist in understanding the evolution of 
the progenitor. The most appealing scenario for K648 is that 
the progenitor was a close binary system that experienced co¬ 
alescence or tidal disruption while both stars were in the MS 
stages and one emerged as a new star with a mass of < 1.6 Mq, 
which then started a new life as the progenitor of K648. This 
progenitor passed through the AGB phase stage, and finally 
became the presently observable C-rich PN K648. 
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Table A 

Detected nebular lines and identification of the HDS spectra. 


•^obs 

(A) 

Ion 

-''lab 

(A) 

Comp. 

/(A) 

/(A) 

5/(A) 

“^obs 

(A) 

Ion 

Alab 

(A) 

Comp. 

/(A) 

/(A) 

5/(A) 

3655.29 

H37 

3656.66 

1 

0.336 

0.043 

0.010 

4879.13 

[Fe III] 

4881.11 

1 

-0.005 

0.051 

0.004 

3655.88 

H36 

3657.27 

1 

0.336 

0.074 

0.013 

4920.01 

He I 

4921.93 

1 

-0.016 

1.290 

0.005 

3656.55 

H35 

3657.92 

1 

0.336 

0.031 

0.011 

4929.28 

[OIII] 

4931.23 

1 

-0.019 

0.038 

0.005 

3657.20 

H34 

3658.64 

1 

0.336 

0.114 

0.015 

4956.85 

[OlII] 

4958.91 

1 

-0.026 

26.803 

4.845 

3657.98 

H33 

3659.42 

1 

0.336 

0.130 

0.016 

4957.06 

[OIII] 

4958.91 

2 

-0.026 

48.170 

5.549 

3658.89 

H32 

3660.28 

1 

0.335 

0.202 

0.018 




Tot. 


74.974 

7.367 

3659.77 

H31 

3661.22 

1 

0.335 

0.190 

0.012 

4970.02 

[OIII] 

4958.91 

1 

-0.029 

0.011 

0.003 

3660.84 

H30 

3662.26 

1 

0.335 

0.275 

0.023 

4975.43 

Ov? 

4977.25 

1 

-0.030 

0.055 

0.003 

3661.97 

H29 

3663.40 

1 

0.335 

0.379 

0.021 

5004.78 

[OIII] 

5006.84 

1 

-0.038 

74.252 

5.601 

3663.26 

H28 

3664.68 

1 

0.334 

0.337 

0.025 

5004.96 

[OIII] 

5006.84 

2 

-0.038 

153.010 

7.626 

3664.66 

H27 

3666.10 

1 

0.334 

0.437 

0.023 




Tot. 


227.263 

9.462 

3666.27 

H26 

3667.68 

1 

0.334 

0.366 

0.023 

5029.97 

Cll 

5032.13 

1 

-0.044 

0.071 

0.005 

3668.03 

H25 

3669.46 

1 

0.334 

0.421 

0.026 

5033.87 

Cll 

5035.94 

1 

-0.045 

0.035 

0.004 

3670.04 

H24 

3671.48 

1 

0.333 

0.486 

0.027 

5045.78 

He I 

5047.74 

1 

-0.048 

0.187 

0.003 

3672.31 

H23 

3673.76 

1 

0.333 

0.535 

0.019 

5059.64 

NIV? 

5061.62 

1 

-0.051 

0.039 

0.003 

3674.93 

H22 

3676.36 

1 

0.332 

0.619 

0.019 

5119.81 

Cll 

5121.83 

1 

-0.065 

0.033 

0.004 

3677.91 

H21 

3679.35 

1 

0.332 

0.697 

0.020 

5515.51 

[Cl III] 

5517.72 

1 

-0.145 

0.021 

0.003 

3681.36 

H20 

3682.81 

1 

0.331 

0.713 

0.020 

5535.63 

[Cl III] 

5537.89 

1 

-0.149 

0.028 

0.003 

3685.38 

H19 

3686.83 

1 

0.330 

0.840 

0.029 

5752.41 

[Nil] 

5754.64 

1 

-0.185 

0.043 

0.002 

3690.11 

H18 

3691.55 

1 

0.329 

1.001 

0.028 

5873.37 

He I 

5875.62 

1 

-0.203 

14.834 

0.138 

3695.70 

H17 

3697.15 

1 

0.328 

1.216 

0.026 

5907.09 

Si I 

5909.37 

1 

-0.208 

0.013 

0.002 

3702.42 

H16 

3703.85 

1 

0.327 

1.407 

0.031 

5977.50 

SII 

5979.76 

1 

-0.218 

0.016 

0.001 

3703.59 

He I 

3705.14 

1 

0.327 

0.722 

0.026 

6148.93 

ClI 

6151.27 

1 

-0.242 

0.045 

0.003 

3710.52 

H15 

3711.97 

1 

0.325 

1.534 

0.027 

6234.73 

Cl 

6237.23 

1 

-0.254 

0.022 

0.001 

3720.48 

H14 

3721.94 

1 

0.323 

1.937 

0.033 

6236.13 

Fell 

6238.39 

1 

-0.254 

0.014 

0.001 

3724.63 

[On] 

3726.03 

1 

0.322 

7.729 

0.156 

6236.39 

Nell 

6238.92 

! 

-0.254 

0.023 

0.001 

3724.61 

[Oil] 

3726.03 

2 

0.322 

9.654 

0.255 

6257.02 

Cll 

6259.56 

1 

-0.257 

0.009 

0.003 




Tot. 


17.383 

0.299 

6257.37 

Cii 

6259.56 

1 

-0.257 

0.013 

0.001 

3727.40 

[Oil] 

3728.81 

1 

0.322 

3.047 

0.166 

6297.90 

[Oi] 

6300.30 

1 

-0.263 

0.225 

0.008 

3727.36 

[Oil] 

3728.81 

2 

0.322 

6.340 

0.309 

6298.44 

[Oi] 

6300.30 

2 

-0.263 

0.041 

0.007 




Tot. 


9.387 

0.350 




Tot. 


0.266 

0.011 

3732.90 

H13 

3734.37 

1 

0.321 

2.428 

0.040 

6309.60 

[S III] 

6313.10 

1 

-0.264 

0.119 

0.005 

3748.69 

H12 

3750.15 

1 

0.317 

2.968 

0.043 

6361.31 

[Oi] 

6363.78 

1 

-0.27! 

0.068 

0.003 

3769.16 

Hll 

3770.63 

1 

0.313 

3.934 

0.052 

6459.20 

Cii 

6462.04 

1 

-0.284 

0.076 

0.007 

3796.41 

HIO 

3797.90 

1 

0.307 

5.291 

0.067 

6459.61 

Cii 

6462.04 

2 

-0.284 

0.022 

0.005 

3818.13 

He I 

3819.60 

1 

0.302 

0.911 

0.015 




Tot. 


0.098 

0.009 

3833.89 

H9 

3835.38 

1 

0.299 

7.105 

0.089 

6522.88 

Nell 

6525.59 

1 

-0.293 

0.014 

0.001 

3867.22 

[Nelll] 

3869.06 

1 

0.291 

9.939 

0.124 

6545.56 

[Nil] 

6548.04 

1 

-0.296 

0.901 

0.012 

3887.32 

H8 

3889.05 

1 

0.286 

22.451 

0.456 

6559.94 

H3 

6562.82 

1 

-0.298 

119.758 

3.906 

3917.39 

Co 

3918.97 

1 

0.279 

0.078 

0.011 

6560.45 

H3 

6562.82 

2 

-0.298 

162.641 

4.167 

3919.12 

Cll 

3920.68 

1 

0.279 

0.167 

0.012 




Tot. 


282.399 

5.712 

3925.02 

He I 

3926.54 

1 

0.277 

0.137 

0.013 

6575.48 

Cii 

6578.05 

1 

-0.300 

0.692 

0.011 

3963.20 

He I 

3964.73 

1 

0.267 

0.780 

0.015 

6580.77 

[Nil] 

6583.46 

1 

-0.300 

1.488 

0.032 

3965.89 

[Ne III] 

3967.79 

1 

0.267 

3.147 

0.043 

6580.93 

[Nil] 

6583.46 

2 

-0.300 

1.692 

0.025 

3968.60 

H7 

3970.07 

1 

0.266 

10.703 

0.167 




Tot. 


3.180 

0.041 

3971.08 

ClI 

3972.45 

1 

0.265 

0.144 

0.015 

6604.63 

Nell 

6607.40 

1 

-0.303 

0.023 

0.003 

3997.69 

Cm 

3999.64 

1 

0.258 

0.035 

0.008 

6628.07 

Oiv 

6630.70 

1 

-0.307 

0.010 

0.001 

4007.72 

He I 

4009.26 

1 

0.256 

0.191 

0.014 

6675.56 

He I 

6678.15 

1 

-0.313 

4.114 

0.055 

4024.64 

He I 

4026.18 

1 

0.251 

1.958 

0.024 

6712.28 

Nil 

6714.99 

1 

-0.318 

0.032 

0.003 

4074.57 

[SII] 

4076.35 

1 

0.237 

0.111 

0.031 

6713.90 

[Sll] 

6716.44 

1 

-0.318 

0.087 

0.002 

4079.83 

Olll 

4081.00 

1 

0.235 

0.034 

0.008 

6724.89 

Cm 

6727.48 

1 

-0.319 

0.037 

0.003 

4100.13 

H6 

4101.73 

1 

0.230 

26.339 

0.248 

6728.30 

[Sii] 

6730.81 

1 

-0.320 

0.133 

0.003 

4119.23 

He I 

4120.81 

1 

0.224 

0.202 

0.013 

6731.54 

He I 

6734.08 

1 

-0.320 

0.022 

0.003 

4142.10 

He I 

4143.76 

1 

0.217 

0.156 

0.006 

6739.57 

Cm 

6742.15 

1 

-0.321 

0.041 

0.004 

4265.50 

ClI 

4267.18 

1 

0.180 

0.660 

0.016 

6741.76 

Cm 

6744.39 

1 

-0.322 

0.061 

0.003 

4265.51 

ClI 

4267.18 

1 

0.180 

0.726 

0.013 

6777.39 

Cii 

6780.60 

1 

-0.326 

0.059 

0.003 

4338.77 

H5 

4340.46 

1 

0.157 

46.674 

0.307 

6798.18 

Cii 

6800.68 

1 

-0.329 

0.052 

0.006 

4361.49 

[Olll] 

4363.21 

1 

0.149 

2.782 

0.026 

6931.26 

He I 

6933.89 

1 

-0.347 

0.053 

0.003 

4386.21 

He I 

4387.93 

1 

0.142 

0.517 

0.010 

7034.62 

Cm 

7037.25 

1 

-0.361 

0.075 

0.003 

4435.82 

He I 

4437.55 

1 

0.126 

0.088 

0.013 

7059.60 

He I 

7062.28 

I 

-0.364 

0.026 

0.003 

4469.76 

He I 

4471.47 

1 

0.115 

4.914 

0.028 

7062.49 

He I 

7065.18 

1 

-0.364 

5.909 

0.100 

4636.98 

Oil 

4638.86 

1 

0.064 

0.055 

0.006 

7092.65 

Sii 

7095.49 

1 

-0.368 

0.015 

0.005 

4640.01 

Oil 

4641.81 

1 

0.063 

0.033 

0.003 

7096.05 

NiV? 

7098.60 

1 

-0.369 

0.020 

0.002 

4701.53 

[Fe III] 

4701.53 

1 

0.045 

0.022 

0.003 

7132.99 

[Ar III] 

7135.70 

1 

-0.374 

0.384 

0.006 

4706.34 

N II 

4708.28 

1 

0.043 

0.030 

0.003 

7278.53 

He I 

7281.39 

1 

-0.393 

0.598 

0.008 

4711.34 

He I 

4713.14 

1 

0.042 

0.672 

0.005 

7316.27 

[Oil] 

7319.14 

I 

-0.398 

0.452 

0.020 

4787.73 

[Eli] 

4789.45 

1 

0.020 

0.110 

0.004 

7317.32 

[On] 

7320.19 

1 

-0.398 

1.340 

0.015 

4859.43 

H4 

4861.33 

1 

0.000 

100.000 

0.188 

7326.89 

[On] 

7329.76 

1 

-0.400 

0.745 

0.015 

4867.28 

[Eli] 

4868.99 

1 

A).002 

0.030 

0.003 

7327.95 

[On] 

7330.82 

1 

-0.400 

0.705 

0.015 







